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Abstract 

We calculate isospin breaking in pion-nucleon scattering in the threshold region in the framework 
of covariant baryon chiral perturbation theory. All effects due to quark mass differences as well 
as real and virtual photons are consistently included. As an application, we discuss the energy 
dependence of the triangle relation that connects elastic scattering on the proton n^p — ► 7r ± p with 
the charge exchange reaction ir~p — > 7r°n. 
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1 Introduction and summary 



In the Standard Model isospin violation (IV) is driven by strong and electromagnetic interactions, 
that is by the differences in the light quark masses rrid — m u and charges Q u — Qd = e, 

£qcd = md 2 mu (uu - dd), £^ ED = ^ («4« - d4d) ■ (1.1) 

Pion-nucleon scattering in the threshold region is particularly well-suited to test our understanding 
of isospin violation because it can be analyzed within chiral perturbation theory (ChPT), which is 
the effective field theory of the Standard Model, and in contrast to pion-pion interactions the quark 
mass difference md — m u already contributes at leading order in the symmetry breaking [1]. This 
problem was addressed in the framework of heavy-baryon chiral perturbation theory in a series of 
papers about a decade ago [2-6]. In [7] we presented a novel analysis of all irN scattering lengths 
based on a covariant formulation of baryon chiral perturbation theory. In this work, the analysis of 
isospin violation in the ttN system will be extended beyond threshold. 

The size of the isospin-breaking effects extracted from low-energy data in early model-based anal- 
yses [8,9] are in conflict with those found in heavy-baryon ChPT [6], while more recently the ChPT 
picture seemed to be confirmed by a tree-level coupled-channel X-matrix approach [10]. The aim 
of this work is to revisit isospin violation in irN scattering in the framework of covariant baryon 
ChPT [11], where the following differences to previous analyses within ChPT should be noted: in [6] 
the full C(p 3 )-amplitude was fitted to data in order to determine the low-energy constants (LECs), 
which are then used to switch off electromagnetic effects to isolate strong isospin breaking. This pro- 
cedure is questionable in the sense that strong forth-order terms are known to be important [12, 13], 
such that electromagnetic LECs might get contaminated accounting for these contributions. For this 
reason, we restrict ourselves to isospin-breaking shifts in the amplitude quantified in terms of the 
so-called triangle relation. Hence, our results will comprise both strong and electromagnetic isospin- 
violating effects. In addition, we include bremsstrahlung (in the soft-photon approximation), such that 
finite terms neglected in [6] are taken into account consistently. Note that here we do not attempt to 
determine the appearing LECs by a fit to the pertinent partial waves, but rather follow the path laid 
out in [7] and collect these from other determinations, concentrating on the isospin-breaking shifts 
induced in the various partial waves. Finally, we hope that the amplitude documented in the present 
work will prove valuable for a consistent incorporation of isospin breaking in forthcoming analyses of 
ttN scattering based on Roy-Steiner equations. 

The main findings and conclusions of the present investigation can be summarized as follows: 

i) In this work, we have systematically analyzed isospin breaking in irN scattering in all mea- 
sured channels (ir^p — ► tt^p, 7r~p — > 7r°n) and the energy dependence of the triangle relation 
in S- and P-wave projections of the scattering amplitude in the framework of covariant baryon 
chiral perturbation theory, accounting for all sources of isospin violation including virtual pho- 
tons and bremsstrahlung. In all cases, we have provided a detailed estimate of the theoretical 
uncertainties. 

ii) Above threshold, isospin violation in terms of the triangle relation amounts to about (2.5 — 
4.0)% between y/s = (1.08 — 1.14) GeV in the S'-wave, the bulk of which is due to virtual 
photons, whereas isospin breaking in the P-waves is very small, at the order of 1 % at most. 
We have compared our work to former studies in heavy-baryon ChPT and analyses based on 
phenomenological models. 
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Our work is organized as follows. Section [2] contains the calculation of the isospin-breaking correc- 
tions to third order in the chiral expansion. After some preliminary chapters on the kinematics and 
isospin structure of the scattering amplitude, we systematically display the calculation of the tree and 
loop graphs, including the effects of virtual photons and soft bremsstrahlung. Many details of these 
calculations are given in the various appendices. In Sect. [3] we analyze isospin violation above threshold 
by concentrating on the triangle relation in the S- and P- waves. It vanishes in the isospin limit and 
is thus an excellent measure of the effects of isospin breaking. We give a thorough discussion of the 
theoretical uncertainty of our calculations. We also compare to earlier work based on heavy-baryon 
ChPT and more model-dependent analyses. 

2 Scattering amplitude 

2.1 Kinematics and isospin structure of ttN scattering 

In this section, we consider various formal aspects of the reaction 

N(p)+AQ)^ N(p')+7r a (q>). (2.1) 

The masses of the initial (final) nucleon and pion are denoted by m\ (mf) and M\ (Mf), respectively; 
m p , m n , M w and M n o are the masses of proton, neutron, and charged and neutral pion, in order. We 
define the isospin limit by the charged particle masses m p and M n . The deviations can be expressed 
in terms of 

A, = M 2 -M 2 , A N = m n -m p . (2.2) 

Furthermore, we will use 

S = p + q = p' + q', A = p — q' = p' — q, A = p — p' = q' — q, Q = p + p' ', (2.3) 

which are related to the usual Mandelstam variables by 

S 2 = s, A 2 = u, A 2 = t. (2.4) 

For convenience, we will often refer to these quantities in the center-of-mass system (CMS) defined by 
p + q = p' + q' = and p • p' = |p| |p'| cos #cms = |p| \p'\z, such that 

_ s + mf - _ s - mf + M 2 _ s + mf - M f 2 _ s - mf + M f 2 

p= %Ts ' q= V5 ' p ' = Vi ' q ' = 2^s ' 



A 1 / 2 (s,mf,M?) , „ , „ A 1 / 2 (s,mf,M?) 

|p| = |q| = — — ' i p " = i q i = — — ' 

t = -(Mf-M^- X ^ M 'K S -^; M ^ -mf + Mf-Mf) 

A 1 / 2 {s,mf,M?) A 1 / 2 ( S ,m f 2 ,M 2 ) m , =mu j h=Ml A(s,m 2 ,M 2 ) _ 

2s Z 2s ^ Z)l 1 ' ' 

where the Kallen function is given by 

A(a, b, c) = a 2 + b 2 + c 2 - 2ab - 26c - 2ac. (2.6) 
u follows from the Mandelstam relation 

s + t + u = mf + mf + M? + M 2 . (2.7) 
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Finally, we will need 

v= j— — v (2.8) 
2 (mi + mf) 

Simplifications of the kinematics in loops will be addressed in Sect. 12.81 In particular, it is convenient 
to introduce a generic nucleon mass m for the 0(p 3 ) contributions, which may be identified with m p 
in the end (cf. Sect. 12 . 5|> . 

As pointed out in [13], the standard decomposition of the pion-nucleon scattering amplitude 

T nN = u(p')^A{s,t) + ^ (t + f)B(s,t)\u(p), u(p')u(p') = 2m h u(p)u(p) = 2m h (2.9) 

is not the appropriate form for a low-energy expansion, since A(s, t) does not fulfil chiral power 
counting due to a cancelation between the leading pieces of A(s,t) and B(s,t). Therefore, T n N is 
parameterized in terms of the two amplitudes D(s,t) and B(s,t) according to 

T wN = u(p'){D(s,t) - l - A$',<j\B{s,t) }u(p), D(s,t) = A(s,t) + uB(s,t). (2.10) 

[ 2(mi + mf) f r J 

The low-energy expansion of D(s,t) and B(s,t) starts at 0(p) and 0(p~ 1 ) and is determined up to 
0(p 3 ) and <D(p) by our leading-loop calculation, respectively. 
In the isospin limit, T n N m ay be decomposed as 

T ab = X ] N , {T + 5 ab + T~\[T\T b ]^XN, (2.11) 

where r a are the Pauli- matrices and xn (xn') is the isospinor for the incoming (outgoing) nucleon. In 
the Condon-Shortley phase convention, the physical amplitudes are related to by (we only display 
the channels to be considered later) 

+ =T + —T~ 

■+p ^ - 1 ! 

-o n = -v / 2T". (2.12) 
2.2 Partial wave expansion of the amplitude 

Above threshold, the partial wave decomposition of the ttN scattering amplitude takes the form [14] 

oo 

T nN {s,t) = V^V^xj { i( l + + lTi-(s)) Pi(z) 

1=0 

-m^-T^ivtfx^PK^Xi, P[(z) = ^P l (z), q=A, (2.13) 

where the Ti± have definite orbital momentum I and total angular momentum j = I ± 1/2. In the 
following, we consider the partial waves fi± defined by 

T l± (s) = ^ i . = 8 J^M, (2.14) 



T - 


= T - 
L 7T 


rp 

± TT+p 


= T n + 


rpCCX 
TT~p 


= T - 

L 7T 
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choosing the normalization such that /o+ corresponds to the S'-wave scattering length. The alternative 
representation of the scattering amplitude 

T wN (s, t) = v / 2^\/2m^A/'iA/'fX f t {g(s, t) + ia • (q x q) h(s, t)} Xi, 



is especially suitable for a calculation in heavy-baryon ChPT. In a non-relativistic expansion, h(s, t) 
and g(s, t) can be interpreted as the spin-flip and spin-non-flip part of the amplitude, respectively. 
The inversion of (|2.13p and (|2.15|) reads 



x/Er, + m\ x /E n i + mr f 

fi±(s)= V P 1{ £ p J ^{g{s^{z) + h{s,t)H\c(\{P l±x {z)-zP l {z))}, (2.16) 

-l 

which can be rewritten in terms of A(s,t) and B(s,t) as 

+ v/ ^T^ , " mf |-A, ±1 ( S ) + f ^±^1 B l±l{ s) \ , (2.17) 



a/E'p + m\JE v , + m f f / mi + m f , _ 
fl±(s) = V TaZ~TT <Ms)+[Vt o B M 



16tt^/s 

where we have defined the partial wave projections 

l l 
Ms) = J dzA(a,t)Pi(z), = J dzB(a,t)Pi(z). (2.18) 

-l -l 

In the isospin limit, (|2,17p coincides with the expressions given in [15]. As pointed out in [5], it is 
favorable not to work with the P- waves P\{s) = fi-(s) and Ps(s) = fi+(s) directly, since in the 
Pi-projection the isospin conserving amplitude itself becomes very small, such that the expected large 
(relative) isospin breaking effect is experimentally not amenable. As we will see in Sect. 13.2.31 the 
results in this projection are very sensitive to the precise values of the LECs and therefore meaningful 
quantitative statements cannot be madeP 1 ! Instead, we use Q(s) and H(s) defined by 

P 3 (s) = g(s)+H(s), P 1 (s) = g(s)-2H(s). (2.19) 

Actually, these projections exactly correspond to the spin-flip and spin- non-flip amplitudes 

i 



G(s) = Y0~^ / dz 9{s,t)Pi{z), 



H{s)= y p w J^_ p 1 J dzh(s,t)\ q \\^\(P 2 (z)-zP 1 (z)). (2.20) 



#1 This is similar to isospin-breaking corrections to the n°N scattering lengths and their normalization with respect to 
the isoscalar scattering length a + , cf. [7]. 
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2.3 Effective Lagrangian and chiral counting 

The framework to systematically analyze isospin breaking in pion-nucleon scattering is baryon chiral 
perturbation theory. In this section, we briefly lay out the framework. Our starting point is the chiral 
effective Lagrangian, 

C cS = + 4 2 i + C® + 4 2 ) + 4 4) + • • • , (2-21) 

formulated in terms of the asymptotically observable pion and nucleon fields, and including also scalar 
sources and virtual photons to incorporate the effects of strong and electromagnetic isospin breaking. 
The superscript in (|2.2ip denotes the chiral dimension, the ellipsis represents terms of higher order not 
needed in our investigation. The couplings appearing in £ e g- are called low-energy constants (LECs), 
they encode information about the higher mass states of QCD that are not active degrees of freedom 
of the effective field theory. 

All monomials in the pertinent fields and sources that make up the terms of the effective Lagrangian 
are constructed according to the following counting rules: 

M n ~<D(p), m~0(l), e~0(p), m u ,m d ~0(p 2 ), t~0(p 2 ), s-s thT ~0(p), (2.22) 

where p is a genuine small parameter and sthr denotes the threshold energy squared of the correspond- 
ing channel. Note that pion four- and nucleon three-momenta are of 0(p). In what follows, we work 
to leading order in isospin breaking, collected in the small parameter 5, 

5 = {m d - m u , e 2 } . (2.23) 

The precise form of the effective Lagrangian employed is given in App. El 

Based on the effective Lagrangian, the pion-nucleon scattering amplitude acquires a chiral expan- 
sion in terms of tree and loop graphs. In this paper, we work to third order in the chiral expansion. 
All the various contributions will be developed in detail in the following sections, here we only make a 
few general remarks. Tree diagrams start at 0{p) (including one-photon-exchange graphs) and obtain 
strong and electromagnetic corrections at orders 0{p 2 ) and 0(p 3 ). One- loop graphs that perturba- 
tively restore unitarity start contributing at 0{p ?J ). Since we work in a covariant formulation of baryon 
chiral perturbation theory, we utilize the infrared regularization method developed in [11] to deal with 
the power counting violations generated by the nucleon mass, m ~ A x with A x = AirF^ ~ 1 GeV the 
scale of chiral symmetry breaking. This method allows for a consistent power counting by uniquely 
separating the soft from the hard parts in any given one-loop diagram. In a nutshell, all one has to do 
is a rearrangement of the Feynman parameter integration corresponding to the combination of light 
and heavy degrees of freedom, 

/ ds(...)-> / dz(...)- dz(...) = I + R, (2.24) 
Jo Jo Jl 

where the infrared singular part I contains the soft (chiral) physics and obeys power counting while 
the regular part, generated from the loop momenta of the order of the nucleon mass, leads to power 
counting violations. However, it can be shown that the contributions to R are polynomials in ex- 
ternal momenta and the pion mass squared. They may therefore be absorbed in the low-energy 
constants of the effective Lagrangian and need not to be considered explicitly. Moreover, the infrared- 
regularization prescription amounts to resumming 1/m insertions in certain heavy-baryon diagrams. 
As a consequence of these higher-order contributions, exact renormalization no longer works: a resid- 
ual dependence on the renormalization scale /U is generated. Finally, the central values of our results 
will be given for a scale \x = 1 GeV, while the sensitivity to the variation of [i can be regarded as 
an indication how large the impact of higher orders will at least be. For more details on the various 
formulations and facets of baryon chiral perturbation theory, we refer to the recent review [16]. 
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Figure 1: Feynman diagrams for ttN scattering at 0(p). Solid, dashed, and wiggly lines, denote 
nucleons, pions, and photons, respectively. 



2.4 Leading-order result (tree graphs) 

In this work, we derive the full amplitude for the channels ir^p — > ir^p and it~p — * 7r°n (referred 
to as the charge exchange reaction T cex = X^fp, which are the only ones directly accessible in ttN 
scattering experiments. 

First of all, we write down the leading-order amplitudes, which will be needed for future reference. 
For this purpose, all nucleon masses may be replaced by m. The diagrams contributing at 0{p) are 
depicted in Fig. [T] and yield within the given simplification 

rLO r , _ _L _ g 2 s + 3m2 _ ^! rLO r , _ _±_ , 9 2 u + 3m 2 2f_ 
a v- P V^) 2p2 2 F 2 s-m 2 t 1 n * + P^' l> 2F 2 ^2F 2 u-m 2 t' 

rLO, t) _ ^2 V2g 2 ( s + 3m 2 u + 3m 2 \ LO _ LO 

B ce AS,t) -~2p2+ \ s - m 2 + u _ m 2 J ' D ceAS,t) - vB cex ( S ,t), 

D L n ° p (s,t) = ^ + uB^ p (s,t), D%(s,t) = ^ + (2.25) 

In the context of bremsstrahlung and the cancelation of infrared divergences we will need the spin 
averaged squared matrix element (T w n = u (p') T n Nu(p)) 

\ M *n\ 2 = \Y, l T ^| 2 = \ Tt {(/ + m (t + m V^ 7 °} » ( 2 - 26 ) 

spins 

for which we find at leading order 

\M%? = (Am 2 - t) (D%) 2 + 2uW%B% -\{t- AM 2 + A, 2 } (b% 
l-M c L e°x| 2 = (4m 2 - t) (-Dcex) 2 + 2i*d£b£ - \ {A (t, M 2 ,M 2 ) + At, 2 } (B^) 2 . (2.27) 

2.5 Tree diagrams beyond leading order 

As depicted in Fig. [2J there are four types of tree diagrams involving vertices beyond leading order: 
contact-term corrections to the Weinberg-Tomozawa (WT) term, to the ttNN coupling and the one 
photon exchange, and mass insertions for the intermediate nucleon. The effect of the latter can be 
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Figure 2: Contact-term contributions. Heavy dots/squares refer to insertions suppressed by one/two 
chiral orders, respectively. 



accounted for by shifting the bare mass m in the Lagrangian according to 

e 2 F 2 

m^m p = m- 4ciM f0 + 2Bc 5 (m d - m u ) — (/i + f 2 + h) 




{2Ml + Ml )+O{p i ) 



{2Ml + Ml )+O{p A ) 



(2.28) 



for the proton and neutron, respectively, where I n is defined in App. IB.2I Actually, we must carefully 
distinguish between the bare mass m in the Lagrangian, which is nothing but the mass of the nucleon in 
the chiral limit, and the masses of the proton and neutron entering the calculation by the requirement 
that the external particles are on their physical mass shell. However, in loops we can still use m as 
a generic nucleon mass to be identified with m p in the end, since the difference m p — m only starts 
at second chiral order and is therefore beyond the accuracy we are working at. For the same reason, 
we may replace m by m p directly in the case of the contact-term contributions. Nevertheless, once 
the shift (|2.28p in the mass parameters of the Lagrangian is performed, the error which is made when 
identifying m — > m p is only of order C(p 3 ). 

\ \ as well as £^ p ^ (see App. [A] for precise definitions), yield corrections to the WT 
vertex. We will always consider s and t as the independent kinematical variables and u to be fixed 
by the Mandelstam relation (12. 7p . Amplitudes originating from WT-type topologies and nucleon-pole 
diagrams will be referred to as vector- and axial-type contributions, respectively. The full vector-type 
tree level amplitudes D v (s,t) and B v (s,t) read 



/7T^cex( S ) ^) 




s,t) 



+ 




2v (di + d 2 ) 
F 2 



2F 2 



v 



32m 3 F 2 i ( 2 ( mp + ra ^ >l> ~ ^ 2 ~ { t + m l- m l) 2 } ( 2 ( m P + m ^) v ~ A tt) 
4vM% d,5 m p + m n e 2 g 7 



F 2 m n 8m 



+ g (2(m p + m n )i/ + m p - m n ) , 



p 



RV / i \ _ _j_ 2m p c 4 4m p t/(d 14 - di 5 ) J_ _L (m p + m n )c 4 

rv /■„ ^ 1 2m p c 4 4m p i/ (d 14 - rf 15 ) 

Corrections to the ttNN coupling are generated by £jj and £^ . They may be incorporated 
replacing g by 

g = g + 2M 2 (2d w - d 18 ) + e 2 F 2 ( 9l + g 2 ) , (2.30) 
g p = g + 2M 2 (2d 16 - d ls ) + e 2 F 2 ( gi + g 2 + g 3 + 9a) - 2B (m d - m u ) (2d 17 - d 18 - 2d 19 ) , 



g n = g + 2M 2 (2d 16 - d 18 ) + e 2 F 2 ( 9l + g 2 - g 3 ) + 2B (m d - m u ) (2d 17 - di 8 - 2d 



19 J i 

for pn, pp and nn couplings, respectively. In terms of these quantities, we find for the full axial-type 
tree level amplitudes A a, (s,t) and B a, (s,t): 

g 2 , ,s-ml Amg 2 s-m 2 2 

A n' P ( S ) = ^2 (™P + mn )j3^ " 2F 2 ( S _ m 2)2 ( S + ^ + 2m ^) > 

g 2 S + mg + 2m p m n Amg 2 a + m p m n 

= ^ +p ( s ,t) = -fljf_ p (u), 

^^cex(s^) = ~ jf^~ s { ( S ~ m n) ( m P + m n) + (s + 2m n m p + 771, 2 ) } 

Amg 2 J s + 2m p m n + m 2 A N m n (3s + m 2 ) + m p (s + 3m 2 ) 



4F 2 1 s-m 2 2 (s-m 2 ) 2 

+ 4F 2 u - m 2 I ^ ~ (mp + mn) ~ if ( U + 2mnm P + m p) } 
Amg 2 f « + 2m p m n + m 2 A N m p (3u + m 2 ) + m n (« + 3m 2 ) 



4F 2 \ u-m 2 2 (u-m 2 ) 2 

1 M , . 99n s + m 2 + 2m p m n Amg 2 s(3m n + m p ) + m 2 (3m p + m n ) 



where 



V2 ccxV ' ' 4F 2 S -m 2 4F 2 (s-m 2 ) 2 

gg p u + m 2 + 2m p m n Amg 2 w(3m p + m n ) + m 2 (3m n + m p ) 
+ 4F^ ~W 2 (u-m 2 ) 2 ' ^ 

Am = m p/n - m p/n = ^ (2M 2 4 (m 2 ) + M 2 J n o (m 2 )) = _^i^±^M + Q(p% (2.32) 

and g 2 , gg p , and gg n , are understood such that only terms linear in the LECs should be retained. We 
have checked explicitly that the double poles in (|2.31|) and the loop diagrams (sx), (si 4 ), and (aio), 
indeed cancel, as required by the analytic properties of the 5-matrix. 
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Finally, \ \ \ and p \ provide contact-term contributions to the one photon 
exchange. The full tree amplitudes D p (s,t) and B p (s,t) turn out to be 

Dl_ p {t) = -Dl +p {t) = — ^- + -j^k + 2e 2 u (de + 2d 7 ) - (h + k 2 + k - 2h 2 + h 4 ) , 

2e 2 2e 2 40e 4 

Bl_ p {t) = -Bl +p {t) = -— (1 + c 6 + or) + j^k ~ (h + k 2 + l 5 - 2h 2 + h 4 ) . (2.33) 

2.6 Renormalization of F and g 

For the numerical evaluation, we aim to replace the pion decay constant F and the axial charge of 
the nucleon g in the chiral limit by the physical values F n and respectively. Since both obtain 
corrections suppressed by two chiral orders, this affects only the leading-order diagrams. The explicit 
expressions in the isospin limit read [17] 

^ = Wl + ^(G-^]og^p)U 0{M% (2.34) 



F 2 V 16vr 2 fa fi 2 
and [13, 18] 

5A - 5 11 + ^ - 16 ^ 2F2 log ^ 2 - 167r2i?2 j + 0(MJ. (2.35) 
Numerically, we use F w = (92.2 ± 0.2) MeV and \g A \ = 1-2695 ± 0.0029 as quoted in [19]. We do not 



consider isospin breaking in F n and <7aCj However, to ensure consistency with the loop contributions, 
the chirally expanded version (|2.35p is replaced by the corresponding expression calculated in infrared 
regularization (IT and I n are defined in App. IB.2j) 

9A = 9 { 1 + [Ml (20m 2 - 13M 2 ) ifto) - 11M 2 /. (m 2 ) - ^ log 

Finally, we remark that the dependence on d\§ cancels between the renormalization of <?a an d the 
corrections to the nucleon-pole diagrams. 



2.7 Wave function renormalization 



As a consequence of the LSZ reduction formula, the full irN scattering amplitude depends on the wave 
function renormalization of pion and nucleon fields according to 



tN 



Z\u{p )T nN u (p) , 



(2.37) 



where T w n is the amputated amplitude, while and Z\ (Z^ and Z^) denote the wave function 
renormalization of the initial (final) nucleon and pion, respectively. 

The Feynman diagrams contributing to the wave function renormalization of pion and nucleon are 
shown in Figs. [3] and HI The corresponding wave function renormalization factors can be deduced from 
the self energy T,(p) via 



Z N = 1 + 



+ 



Zjr 



1 + 



dE(p) 



dp 2 



+ 



p 2 =M% 



(2.38) 



# 2 In fact, radiative corrections are taken into account when extracting F n and g\ from experiment. Moreover, the 
final accuracy of our results will be rather limited by unknown LECs than by the uncertainties in <?a and F„. 
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Figure 3: Wave function renormalization of the charged and the neutral pion. 

/'""*'* \ jJ^^V /''^"\ 



Figure 4: Wave function renormalization of the nucleon. 
where the ellipsis indicates higher orders. Note that the wave function renormalization factors depend 



on the parametrization of the pion fields In the a parametrization (|A.5|) we obtain^ 4 



p 2 M 2 C 9D 

^ = 1-^(^ + 1 + ^)-^- 2 {e*(k 1 + k 2 ) 

M 2 C 7T o 20 , o , 

Z n o = 1 - j^Tff ~ y e + M " 2e 2 (fc 4 - 2A: 3 ) , 

,2 



2 z' 

Z n = 1 + Aj|2 (4m 2 - 2M 2 ) M 2 il(0) + (4m 2 - 2M 2 ,) M 2 ^°(0) - 4M 2 / W (m 2 ) 

M 2 M 2 1 e 2 

- 2M 2 I 7r o (m 2 ) - ^fC, - Y^C.ro h Z p = Z n -^(L 7T + l + C 7T ). (2.39) 



Thus, the wave function renormalization factors are UV and (in the case of the charged particles) IR 
divergent and, with respect to the cancelation of divergences, constitute therefore an indispensable 
part of the amplitude. Their contributions D wi (s,t) and B wi (s,t) due to 5Z = Z — 1 expressed in 
terms of the leading-order amplitudes (|2.25|) read 

Df ±p {s, t) = {5Z V + 5Z P ) D%{s, t), B£ p (s, t) = (5Z n + 5Z P ) B%(s, t), 

D^ x (s, t) = 1 (5Z n + 5Z w o + 5Z P + 5Z n ) D™{s, t), 

Bcl(s, t) = ± {5Z„ + 5Z T o + 5Z P + 5Z n ) B™{s, t). (2.40) 
2.8 Loop diagrams 

At order 0(p 3 ) loop graphs solely contain vertices from the leading-order Lagrangians (|A.2p . The list of 
diagrams depicted in Figs. [SHU] comprises contributions from purely strong diagrams, virtual photon 
corrections to vector-type and axial-type diagrams, and, for the elastic channels, loop corrections 
to the one photon exchange. At the order we are working all nucleon masses may be replaced by 
a generic nucleon mass m (cf. Sect. 12. 5p . which simplifies the calculation tremendously. However, 
counting e ~ 0(p) the pion mass difference is of the same chiral order as the pion masses themselves 
and therefore has to be kept. Apart from virtual photons, this is the second source through which 
isospin violation manifests itself in the loop contributions. 



* 3 At C(p 3 ), is accidentally independent of the parametrization, since no vertices involving more than two pions 
occur. 

* 4 Cr, L n , and the loop functions I a, I it, are defined in Apps. [BTT1 and [fOl respectively. 
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The evaluation of the diagrams proceeds as usual: at first, the Dirac structure is simplified by 
application of the Dirac equation, such that the power of the loop momentum in the numerator 
is reduced as far as possible. The remaining tensor integrals are reduced to scalar loop functions by 
means of Lorentz invariance. As for the strong diagrams, we follow the conventions in [13] for the basis 
vectors to facilitate the comparison to this work. The details of the tensor decomposition can be found 
in App. 1531 where in particular the question of how to properly incorporate the pion mass difference 
is addressed. Since u(p')$'u(p) = u(p')du(p) up to higher orders, A(s,t) and B(s,t) as defined in (|2.9j) 
can be simply read off as prefactors of 1 and ^, respectively, and D(s,t) then follows from (|2.10p . 
The results for the individual diagrams are given in App. O To arrive at these expressions, we have 
made use of the relations between the basis functions, which in some cases simplifies the formulae 
considerably. 

The results for the scalar loop functions calculated in infrared regularization are presented in 
App. IB.2I As most of the work to arrive at these expressions is rather straightforward but tedious, 
we will focus in the remaining part of this section on the issues which are more intricate, namely the 
calculation of the so-called box graphs, the extraction of infrared divergences and the separation of 
the Coulomb pole generated by virtual photons in a threshold expansion of the amplitude. 



2.8.1 Box graphs 

We will refer to any diagram whose loop contains four propagators as "box graph" . These topologies 
lead to scalar loop integrals which are technically the most difficult ones we have to deal with, as 
three Feynman parameters are needed to perform the momentum integration in the usual manner. 
Having done the Feynman parameter integration corresponding to infrared regularization, we are still 
left with two Feynman parameter integrals. The analytic properties of these integrals are far from 
obvious, which in particular causes numerical integration to fail as soon as one goes above threshold. 
Let us first consider the box diagram consisting of three nucleon and one meson propagator 



In(s,t) 



d d k 



i j\ (2nY ! (iUf - /,--') i/n 2 - (p - /.- ) 2 1 ( in' 2 - (S- I:) 2 )(m 2 - (// - A- ) - ) 



(2.41) 



For simplicity, we will neglect the pion mass difference in this section (the general case is treated in 
App. IB.2T5]) . Below threshold, the representation 



Ii3(s,t) = dx dy 



Co + z { 



16tt 2 A 2 I 2c 



3/2 



arccos 



^0 



A = (1 — y)s + ym 2 — x(l 
Ml + (1 - y)(s - m 2 ) 



2A 



x)y 2 t - y(i 

Ml 
A 



y/°o + zl 

y)Ml 



+ 



ZQ 

2c 



Zq, 



(2.42) 



may be used. One way to obtain a representation valid above threshold is the application of dispersion 
relations, since the regular part is numerically harmless. However, Iis(s, t) is also needed with the pion 
replaced by a photon, and the virtual photon leads to an IR divergence. The prescription how to deal 
with the vanishing photon mass is to send m 7 — > at the very end of the calculation. Therefore, we 
first have to analytically continue above threshold and thereafter put m 7 — > 0, which cannot be done 
on the basis of (|2.42p . As shown explicitly in the appendix, dispersion relations are not appropriate 
either, such that it is inevitable to improve (|2.42p . i.e. the y integration should be performed. The 
second term is a rational function in y and can be integrated by means of a partial fraction expansion. 
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Figure 6: Strong one-loop diagrams for tt + p — > 7r + p. 
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Figure 7: Strong one-loop diagrams for tt p — ► 7r°n 
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Figure 8: Electromagnetic vector-type diagrams for it p — > it p. 
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Figure 9: Electromagnetic vector-type diagrams for ir + p — ► 7r + p. 
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Figure 10: Electromagnetic vector-type diagrams for 7r~p — ► 7r°n. 
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Figure 11: Electromagnetic axial- type diagrams for ir p — > it p. 





Figure 12: Electromagnetic axial-type diagrams for n + p — > 7r + p. 
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Figure 13: Electromagnetic axial-type diagrams for it p — > 7r n. 




Figure 14: Photon-exchange diagrams for — > ^p at C(p 3 ). 



19 



After an integration by parts 



dy 



y ( c q + z o) 

2A 2 Cn /2 



arccos 
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4M 2 { 2 y(2/i + 2/2) - 4yiy 2 } arccos 



go 



4Mj 



dy 



2y(yi + 2/2) - %i?/2 



(c 2 ) 3/2 (yi - y2) 2 y / (yi - y)(y - 2/2) 



(c 2 ) 3 / 2 7 (yi - y 2 )V(2/i -2/)(2/-2/2) d 2/ 



arccos 



■20 



(2.43) 



where the notation of App. IB.2.51 has been used, also the first term is reduced to an integral over a 
rational function in y, since as far as square roots are concerned 



dy 



arccos 



z 



\J{y\ -y)(y- 2/2) 



(2.44) 



The x integration turns out to be harmless, such that the continuation above threshold can be per- 
formed and the IR divergence for m 7 — > extracted. The final results are again summarized in 
App.[BJL5j 



2.8.2 Extraction of infrared divergences 

The following scalar loop functions possess an IR divergence: l\(t) (cf. diagram (^3)), V2o(t) (1*4), 
V\i(s) (vi), l! 3 (s,t) (a 7 ), A 2 i(s,t) (ag), Ai 2 (s,t) (a 6 ), and Pn(s,t) (p 9 ). Thereof, the first two 
integrals do not induce any problems and Ii 3 (s,t) has already been treated in the previous section, 
but unfortunately this procedure cannot be directly adopted to the remaining cases. Let us first 
consider Vii(s), which is the simplest remaining task, as A2i(s,t), A\2(s,t), as well as Pn(s,t), are 
box graphs. The following procedure is inspired by [20]. 

We extract the IR divergence from the loop integral calculated in ordinary dimensional regulariza- 



tior * 5 



dxi dx2 dx3 5(1 — x\ — x 2 — X3) 



167T 2 J x\m^ + x\Ml + x\m 2 — X2X3 (s 



The substitution 



X 2 -> (1 - X 1 )x 2 , 



•1'3 



'1 - xi)x 3 



implies 



which leads to 



5 (1 - Xl - X 2 - X 3 ) -> 5 ((1 - Xl)(l - X2 - x 3 )) 



5 (1 - x 2 - x 3 ) 

1 — Xl 



(2.45) 
(2.46) 
(2.47) 



1 



16vr 2 



1 

dx / dy 



1-2/ 



ym 2 + (1 - y) 2 s(x) 



s(x) = xm 2 + (1 - x)Ml - x(l - x)s. (2.48) 





Performing the y integration, we finally arrive at 



1 



32vr 2 



L 







dx 

s(x) 



+ 



dx log 



s{x) 







s(x) 



(2.49) 



* 5 The regular part is finite in the limit of vanishing photon mass, such that m 7 may be put equal to zero right from 
the start. 
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For the remaining box graphs, this approach has to be generalized, which we will illustrate for 
^4i2(s,i) (^2i(s,t) and Pn(s,t) are treated analogously). The starting point is again 

.dim/ + n _ _J_ f dgi dx 2 dx 3 dx 4 5(1 -x 1 -x 2 -x 3 - s 4 ) 
mm^ — X3 (s — m 2 ) + (x 2 </ + + x 4 p) 

Separating (and performing) the x\ integration with the help of the substitution 

Xi ->• (1 - xx)xi, ie {2,3,4}, 5(l-x 1 -x 2 -X3-x±)->— — ^ — — — — (2.51) 

1 — X\ 



yields 



'ill"/',, M— ' f .1 .. .1 .. .1 .. / "'" '' •'' ^ 



^iT(M) = T7TT / dx 2 dx 3 dx 4 



2 — / '° — m 2 



167T 2 J [ {/(^2, %3, Xa) ~X 3 (s- TO 2 )} S 2 (x 2 , X 3 , X A ) 



.of ml+xzis— m 2 \ m 2 .+X3(s-m 2 )—2f(x2,X3,X4) 

4mr, arctan ,' v ' — arctan — '- v , 

' V ^-S 2 (x 2 ,X 3 ,X4) yJ-S 2 (x2,X 3 ,X A ) , , 

(-s^(x 2 ,x 3 ,x 4 ))' ;1 /^ J 
s 2 (x 2 ,x 3 ,x 4 ) = m^ + x 2 3 (s - to 2 ) 2 - 2to 2 {2f(x 2 , x 3 , x 4 ) - x 3 (s - m 2 )} , 

/(aj 2 , x 3 , £4) = (aW + x 3 E + x 4 pf . (2.52) 
The structure of this formula suggests to repeat the same steps for x 3 

x . _» (1 _ X3 ) Xi , i G {2,4}, 6{l-x 2 -x 3 -x 4 ) -> ^-^-^ ^ (2 53) 

1 — x 3 

In this way, the X3 integration factorizes and the 5-function can be applied. The result of this manip- 
ulation remarkably resembles (I2.42h . and indeed the X3 integration can be treated in the same way as 
the y integration in the previous section. Thereafter, the limit to 7 — ► may finally be taken. In this 
way, we find the IR divergence ~ logm 7 and the correct finite terms. As compared to (|B.24|) . putting 
the photon mass to zero at the end simplifies the result considerably, while the derivation requires 
about the same amount of work. 



2.8.3 Coulomb pole 

The ttN scattering amplitude at threshold contains Coulomb singularities once virtual photons are 
included. In particular, the threshold expansion of loop diagrams which allow for an exchange of a 
virtual photon corresponding to a static Coulomb potential involves a term oc l/|p|, where p is the 
CMS momentum of the incoming particles. The coefficient of the Coulomb pole can be associated with 
the ttN scattering amplitude itself. In this section, we gather the parts of the amplitude contributing 
to the Coulomb pole in order to confirm this statement. 

The diagrams which give rise to the Coulomb pole are (v\) and (05) for ir~p — > tt~p, (v±) and (a^) 
for ir + p — > ir + p, and (v\), (05), and (ag), for it~p — ► ir°n, respectively. Based on the representations 
given in App. lB~2.9l and fB.2.1H one can show that for p — > 

32 (to + M n ) |p| 32(TO + M 7r )|p| 

=> lim {|p| (V n (s) + (u - m 2 ) A 12 (u,t))} = 0. (2.54) 
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For Ai2(u,t), this relies on the observation that 

gt P (u, 0) = glTiu, -A.) = ~ (m+ ^ )|p| + 0(1). (2-55) 

We denote the analogues of the tree amplitudes (|2.25p with the one-photon-exchange contributions 
oc 1/t subtracted by L> LO (s,t) and B h °(s,t). The Coulomb pole in Vn(s) may be written as 

<>>- 16 AV Ms 1 ,„ 2 . Ag) = 32( m+ 1 M,)|p| +0 ^')- P' 56 > 

Expressing the pole in A±2(u,t) by Vjj° le (s) via (j2.54[) . we find that the Coulomb-pole contributions 
to the full amplitude can be summarized in the simple form 

Bf ole (s, t) = -2(Q in + Q ont )i e 2 {s - m 2 - M 2 ) V» ole (s)Bf°(s, t), i E {vr^, vr+p, cex} , 
Z?f ole ( S) t) = -2(Q in + Q ont ) ie 2 (s-m 2 - M 2 ) V?° lc (s) b\°{s, t), (2.57) 

where Qm/out are given as the products of the charges of the particles emitting and absorbing the 
photon in the initial and final state, respectively, such that 

(Qin + Qout)ir-p = — 2, (Qin + Qout)7r+p = 2, (Qin + Qout)cex = 1- (2.58) 

It is obvious from the above that the perturbative treatment of photon exchange breaks down for 
very small |p|. The Coulomb pole cx l/\p\ discussed above is the leading, 0(e 2 ) approximation to the 
Gamow-Sommerfeld factor [21,22] 

\ 27T V -, , /nt 2\ Qin/oute 2 S - T71 2 - M 2 

G(T?) = exp(2vr,)-l =1 - 7rT? + 0(??) ' ^ = ^ X l / 2 {s,m 2 ,M^y ^ 

which resums Coulomb photon exchange to all orders. Obviously, the one-photon approximation is 
insufficient as soon as rj ~ 1. If one attempts to describe physical data also very close to threshold, one 
has to multiply the amplitude with a factor of \J G(r/i n )G(ri ont ); in order to avoid double counting, we 
have to remove the approximation to the Coulomb pole as contained in our one-loop representation. 
In Sect. El we will discuss results for amplitudes purified of Coulomb photon effects in this sense. 



2.9 Bremsstrahlung 

The appearance of infrared divergences in T n pj renders the amplitude ill-defined if the photon mass is 
sent to zero. The solution to this problem is well known [23-25]: it is inconsistent to merely account 
for virtual photons, since real photons can be emitted with arbitrarily low energies in the final state 
as soon as charged particles are involved, which inevitably requires the inclusion of bremsstrahlung. 
One can never specify the number of escaped photons exactly, but has to decoherently sum over all 
processes integrating the photon energies up to some typical detector resolution E max (and over the 
phase space available). The integration over propagators from which a photon is radiated generates 
another infrared divergence which cancels the singularities caused by virtual photons. In this way, the 
inclusive cross section becomes infrared finite. 

To the order we are working, it suffices to consider the emission of a single photon to achieve this 
cancelation of infrared divergences in the total cross section 

<7 tot (s,t,£; max ) = a nN (s,t) + a 7 (s,t,E max ). (2.60) 
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Figure 15: Bremsstrahlung diagrams for tt p — > ir p. 



Denoting the photon momentum by k = (E^, k) and the amplitude for bremsstrahlung by 

T 7 = u{p%(k)T^p, P ',q,q',k)u(p), (2.61) 

the cross section for the pure scattering process a n N(s,t) and for the bremsstrahlung correction 
er,y(s, t, E m3X ) can be obtained from the spin and polarization averaged squared matrix elements 

spins, pol. spins 



via 



*M*,t) = J j ^2\M^ N \ 2 {2iTf5^{p + q-p' -q'), 
a 7 (s,t,E max ) = I J \M 7 \ 2 (2n)^(p + q-p' - q' - k), (2.63) 



where the flux factor T and the two-body phase space d</>2 are given by 

^ = V(p-,)»-mV. dfe - ^L-^- . (2.64) 

Since the flux is determined solely by the incoming particles, J- is the same for both processes. The 
Feynman diagrams for bremsstrahlung are depicted in Figs. [T51I17I 



2.9.1 Soft-photon approximation 

It is clear from the previous discussion that in general the infrared divergences only cancel at the level 
of cross sections. In particular, one has to perform the phase space integration of the amplitude, which 
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Figure 16: Bremsstrahlung diagrams for ir + p — > 7r + p. 




Figure 17: Bremsstrahlung diagrams for 7r p — > tt n. 
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is a considerable task if analytical expressions are aimed for. The amount of work can be reduced 
significantly be means of the so-called soft-photon approximation, in which already the sum of the 
squared amplitudes is infrared finite. It neglects terms of order 0(E max ), but reproduces singular 
structures and constant terms correctlyl^ 6 ! 

The crucial approximation is to drop the photon momentum in the delta function. In principle, 
this amounts to reducing the three-body phase space integration to a two-body one by setting k = 0. 
In fact, every term in |M 7 | 2 , which is not of order 0(k~ 2 ) may be neglected, such that effectively only 
expressions of the form 

(2pi ■ k ± m 2 ) (2p 2 ■ k ± m 2 ) ^ ^ 

need to be integrated over the photon momentum. That the error is indeed of order 0(E max ) follows 
from dimensional analysis 

d3k 1 r * k log(*=Y (2.66) 



(2vr) 3 2£ 7 2p 1 -k2p 2 -k J k 3 \ m 7 

d 3 k 1 f d 3 k 



(27r) 3 2E 7 2 Pl ■ k J k 2 

Ery ^i/max 

In particular, it is clear that terms of order 0(1) can only originate from the phase space integration 
of (|2.65p . Moreover, the m 2 terms in the denominator of (|2.65p may be omitted, while the full energy- 
momentum relation E 2 = k 2 + m 2 is necessary to arrive at the correct constant terms. Neglecting rrt 7 
therein merely reproduces the infrared divergence. This result can be understood from the observation 
that k enters the energy-momentum relation quadratically, such that both terms are of the same size 
in the integration region where |k| ~ m 7 , while the m 2 terms in (12.650 are still suppressed by one 
order of m 7 in this domain and therefore constitute a higher-order effect. For p\= P2=V = (E p ,p), 
we obtain 

d 3 k 1 1 ( ml E p E p + \p 



(2vr) 3 2£ 7 (2p ■ k) 2 32vrV 

m ax 

while the formula for the general case reads 



log -m~ + e log + ° {m ^ (2 - 67) 

4 ^max iPl h P ~ iPl J 



J (2vr)32£ 7 2p x ■ k 2p 2 • k 32vr 2 \ A ° g 4^ ax J p 2 + J p 2 |p| i0g E § - |p| J + U[ " h 

&y<-Emax 

p = xp x + (1 - x)p 2 = {E^ p) . (2.68) 

Finally, the soft-photon approximation allows for the application of the usual expressions for the 
Mandelstam variables, since the difference (p 1 + q') 2 — s = (p + q — k) 2 — s = —2k ■ (p + q) + m 2 — > 
—2k • (p + q) yields only contributions of order 0(E m3iX ) even when multiplied by 0(k~ 2 ) terms, which 
is beyond the accuracy of this approximation. 



i l 



#6 A 

comparison of the soft-photon approximation and the exact result in the context of radiative corrections to kaon 
decays can be found in [26]. 
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2.9.2 Bremsstrahlung amplitude 

The contributions from the individual diagrams are given in App. IC.51 Their sum can be brought into 
the following form 

K- P = a 4^i + 4i + 44 + - tf'h" + ^ + 4, H^ p = r% ( q «- - q ') , 

elf 2 4m 2 g 2 \ e if 2 Am 2 g 2 



2F2 2p'-A;[ y (p + g ) 2 -m 2 J 2F 2 2f/:[ y (j/ + q') 2 - m 2 

2F 2 \ 5 (p + g) 2 -m 2 J\ 2q' ■ k 2p>-ky 

e f 2 4mV 1 f (2q - kY 2q» 

l ~9 -j-r—iTo — —o}< K—r. + 



2F 2 { y + g') 2 -m 2 } { 2q ■ k 2p ■ k 

meg 2 ( 1 1 1 e(l - g 2 ) ( 2p ■ q + q 2 2p' • q' + </ 2 



"' " F 2 \2pMfc ~ 2p-fcJ' ° 5 " 2F 2 \ 2p' • k 2p ■ k 

meg 2 f (2</ + jfe)" (2g - fc)^ 2(g - g') A 

F 2 \ 2q'-k 2q ■ k 2p' • k 
\[2meg 2 



°6 



F 2 ({p' — q) 2 — m 2 ){{p — q') 2 — m 2 ) ' 

2J2 



v/2e 1 f 9 2m V 2m 2 « 



2F 2 2p ■ k \ ' (p' + g') 2 — m 2 {p' — q) 2 — m 2 

V2e{2q-kY \ 2 2m 2 g 2 2m 2 g 2 

3 2F 2 2g • k \ 9 (p 1 + q') 2 -m 2 (p - q') 2 - m 2 

V2e 1 |\ 2wo , , , /2 . 2m 2 g 2 



2F 2 2p 



^{(l-^)(2p'.,' + ^)- (?/ _ 2 7_ m2 2,.,-} 



(2.69) 



where terms generated by the one-photon-exchange diagrams are omitted. On the level of the soft- 
photon approximation they can be restored by the replacement 

Ac 2 VP 2 

l-g 2 ^l-g 2 -^f- (2.70) 

in ai, a^, ag and 05. As a check on (12.69p . we have verified explicitly that the Ward identity 

k^_ p = k^ +p = k,T? ex = (2.71) 

is fulfilled. Besides m 2 terms in the denominators, (|2.69p is still exact. The soft-photon approximation 
implies that all /c's in the numerators may be dropped and the usual definitions for the Mandelstam 
variables applied, since we are only interested in 0{k~ 2 ) terms in |A^ 7 | 2 in the end. In particular, this 
means that 61 does not contribute. In view of this procedure, one might worry what happens to gauge 
invariance. Indeed, within the above approximation (|2.7ip is only valid up to terms of order 0(k), as 
the terms omitted in Tf* are by construction of order 0(1). However, using the photon polarization 
sum in the form [27] 

^2 e* fl (k,X)e v (k,X) = -g lJU , + k li c u + k 1/ c fi , c = -7^, c° = — !— , (2.72) 

A=±l ' ' 



26 



it is clear that in 



\ M l\ 2 = \ E ^(fc,A)e,(A:,A)Tr((/ + m)r^^ + m)7 (^)t 7 0) 



(2.73) 



A=±l 



the additional terms due to k^T^ are at least of order 0(k~ 1 ), which goes beyond the soft-photon 
approximation. In this sense, gauge invariance is only maintained up to the order one is working at. 
Nevertheless, it follows from (12 . 721) and (12.730 that the squared matrix element for the bremsstrahlung 
process is still given by 



1 



TV ((/ + m)T»(f + m) 7 °rt 7 ) . 



(2.74) 



The evaluation of the traces in (I2.74|) is now straightforward, but tedious. Since we are eventually 
interested in 

■\M y \\ (2.75) 



\M* 



(27r) 3 2£ 7 

the spin- and polarization-averaged squared matrix element integrated over the photon energies, all 
that remains to be done is the k integration, which may be performed on the basis of (|2.68p . The final 
result can be expressed as a correction factor to the leading-order formula for ttN scattering (|2.27p 

„2 



\M,\ 2 



\M^\ 2 C(s,t,E n 



where 



= log 

+ 2 



mz 



8vr 2 
2m 2 



(2.76) 



AE 2 
s -m 2 - Ml 



-f(t,m) 



2Mi 



t 



gu (s,m,M n ) 



u — m 



Ml 
2 -M 2 



f(t,M n ) 



-gu (u,m,M n ) 



+ 2m 2 g pp (s) + 2Mlg qq {s) - 2(s - m 2 - Ml)g pq {s) 
- (2m 2 - t)g pp ,(s,t) - (2Ml - t)g qq ,{s,t) - 2{u - m 2 - Ml)g pq ,(s,t), 



log 



AE 2 



m 



4 - 
Ml 



2m 2 



t 



f(t,m) 



2Mi 



t 



g u (s,m,M n ) + 2 



U — 777 



M 2 
2 M 2 



f(t,M n ) 



-gu (u,m,M n ) 



" 9 ail V"' ■"■'t/ 1 9 

+ 2m 2 g pp {s) + 2Mlg qq {s) + 2{s - m 2 - Ml)g pq {s) 
- (2m 2 - t)g pp ,(s,t) - (2Ml - t)g qq ,{s,t) + 2(u - m 2 - Ml)g pq ,{s,t), 



C cex (s,t,E max ) = log 



777r 



AE 2 



2 + 



777 



777^ 



M 2 

—^511 (s,m,M v ) 



+ m 2 g pp {s) + Mlg qq {s) - (s - m 2 - Ml)g m {s). 



(2.77) 



f(t,m) and gu (s, m, M n ) are defined in App. IB. 11 while the remaining functions are relegated to 
App.[BlJ Note that the pion mass difference does not play a role for either of the three channels. For 
the charge exchange reaction this is true, as q' does not enter. Furthermore, C ce x (s, t, E max ) is in fact 
independent of t. Therefore threshold kinematics are given by Stm- = ( m + M^) 2 , u tnj - = (m — M n ) 2 
and t = 0, since s is the only Mandelstam variable needed for the charge exchange reaction. Using 

777 . . 777 , _ „ 

511 Khr) = — » ( 2 - 78 ) 



/(0,777) = /(0,M w ) = l, 5l l( Sthr ) 

2 



9pp{ s thr) — 9pp' ( s thrj 0) 



771/ 



9qq(s t hr) = 9qq< (sthri 0) 



Myr 



Mr 



9pq( s thr) — 9pq' ( s thr j 0) 
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one can easily show that C(sthn 0, -EWx) vanishes as it should, since at threshold there is no phase 
space available for the emission of bremsstrahlung any more. 

2.10 Total amplitude 

Eventually, we now have everything at hand to write down the full amplitude and show how the 
cancelation of UV and IR divergences works. First, the results of Sects. [23| 12.71 and l2.8} are combined 
to 

Bi(M) = ^(M) + B?(s,t) +Bf op (s,t), i E {^-p,^+p,cex} , 
A(«,t) = Df°*(a,t) +Df(s,t) + D l ° op (s,t), 
Bt e (s,t) = BJ(s,t) + B?(s,t) + B?(s,t) = Bf (s,t) + O(l), 

Df ee (s,t) = Dj(s,t) + Dl(s,t) + Df(s,t) = D\°(s,t) + 0(p 2 ). (2.79) 
Second, the generalization of (12. 27ft beyond leading order reads 

\M^ p \ 2 = (4m 2 ~ «) D%D% + ut [d%B% + D%B%} -\{t- 4M$ + V} B%B% 

+ 2 (4m 2 - t) D% + D%) + 2vt [ D % (V° ± °J + B%) + B% + D%) } 

-i{t- AMI + 4- 2 } B% (b^J + B%) , 

= ( K + m n) 2 - *) M o + _ A^) {M o + , M 

- ^{a (i, Ml, M 2 ) + 4^ 2 + (2M 2 + 2M ff 2 „ - t) A 2 , - 4 I /A N A 7r }^ e 5 c L ° 

+ 2 (4m 2 - t) D C L ° (z£?> + + 2«/f {d c L ° (s*° p + 5 c w c f x ) + B^ (l& + £>*) } 

- 1{A (t, M 2 , M 2 ) + 4^ 2 }l? c L ° (bJS? + B£) . (2.80) 

To arrive at these expressions, we have used the following: the amplitude has been determined up to 
0(p 3 ). Hence, the resulting squared amplitude can be described consistently up to 0(p 4 ) only, since, 
by interference with the leading order, 0(p 4 ) terms of the amplitude start contributing at 0(p 5 ). 
Therefore, terms of the form ^ loo P^ wf ) etc., are of higher order in the chiral expansion. For the same 
reason, the imaginary parts of D and B were omitted. Strictly speaking, the squared amplitude looks 
like 

\M\ 2 = (4m 2 - t)\D\ 2 + vt (BD* + B*D) - - {t - 4M 2 + 4^ 2 } \B\ 2 , 
but the chiral order of the terms involving imaginary parts is at least 0(p^). 



2.10.1 Cancelation of ultraviolet divergences 

We can only expect the UV divergences to cancel at 0(p 3 ), higher-order singularities as introduced by 
infrared re gula rization have to be removed by hand. Collecting the prefactors of the UV poles of B(s, t) 



and D(s, t wA in B x (s, t) and D x (s, t), the cancelation should work up to 0(p) and 0(p 3 ), respectively. 
The chiral expansion is performed in terms of the rules given in (|2.22p and M 2 = M 2 — 2e 2 F 2 Z. In 
this way, we have checked explicitly that indeed 

B x ~ 0(p 2 ), D x ~ 0(p 4 ). (2.81) 



#7 



Note that we cannot use A(s, t) for this purpose, as it does not fulfil chiral power counting. 
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2.10.2 Cancelation of infrared divergences 



In contrast to the UV divergences, the IR singularities vanish exactly once bremsstrahlung is included. 



Inserting the IR divergent parts of Z p , Z w , and the pertinent loop functions, into (|2.80p . we obtain^ 8 



■>-^- M ~ •,,,„,M,)-2"- m2 - M - 2 

ml ( . 2m 2 -t ., , 2M 2 - t 



+ 2 2 "flu ( s > m > M -) " 2 2 m ' M -) f + 



s-m 2 - Ml , N u-m 2 - Ml , A 

-gn (s, m, M n ) + 2 = ^g n (u, m, M v ) } + 



g2 772,2 /• 2 A/f 2 ^ 

|A^cex| 2 = - i -2l^| 2 log-^h + 5 ^g u (s,m,M n ) +•••, (2.82) 

which actually cancels the logm 7 terms in (12. 76ft . Thus, we have shown explicitly that 

\M^ N \ 2 + \M^\ 2 (2.83) 

is infrared finite. 



2.10.3 Checks on the amplitude 

First of all, the presence of UV and IR divergences in the theory provides a powerful consistency 
check on the amplitude: in the final result both types of singularities must cancel. Beyond that, we 
have shown that the coefficient of the Coulomb pole is indeed related to the ttN scattering amplitude 
(Sect. [278.3P and that at threshold the IR divergences cancel among themselves, while bremsstrahlung 
is absent (Sects. 12X21 and 12.10.21) . 

In order to check the representations of the scalar loop functions based on Feynman parameter 
integrals we have also invoked dispersion relations (cf. Apps. IBTTI and |B~2|) . If the dispersive represen- 
tation does not converge due to an IR divergence, we have at least checked numerically that dispersion 
relations for finite m 7 reproduce the Feynman representation with explicitly separated IR divergence 
when the photon mass is made sufficiently small. 

Finally, we have verified that in the isospin limit the result quoted in [13] is reproduced. 



3 Isospin violation above threshold 

First of all, one should note that the concept of a partial wave expansion becomes doubtful once 
electromagnetic interactions are included, since a partial wave projection of the 1/t terms generated by 
the one-photon-exchange topologies is not feasible. However, as isospin violation due to static Coulomb 
interaction is in principle well understood, we will focus on the real part of the scattering amplitude 
with the one-photon-exchange contributions amputated and the Coulomb pole as separated in (|2.57|) 
subtracted!^ 9 ! This procedure ensures that T n N is finite at threshold and guarantees consistency with 



the treatment of the scattering lengths in [7]. 

* 8 The ellipsis denotes terms finite in the limit m 7 — > 0. Note that due to (|2.78|) the coefficients of logm 7 vanish at 
threshold. 

* 9 To avoid an unnecessarily clumsy notation, the resulting amplitude is again denoted by T^jv. 
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3.1 Triangle relation above threshold 

We will quantify isospin violation above threshold in terms of the triangle relation 

%w .,rt^H!CT. ( 3.D 

/,' ± p (»)-/r± p «+^/rM 

However, this expression involves several problems: first, the IR divergences due to virtual photons 
no longer cancel among themselves, such that the inclusion of bremsstrahlung is inevitable. Never- 
theless, even in the soft-photon approximation this is in principle only possible at the level of squared 
amplitudes. Second, elastic and inelastic channels have a different threshold behavior in the P-wave 
projections, which renders (|3.ip useless near threshold. 

To remedy the bremsstrahlung problem, one might try to define an analogous triangle relation for 
the squared amplitudes, which may be regarded as infrared safe objects. Obviously, there is no way 
to combine 

| T ^+ P |2 = | T + _ T -| 2 ; | T ^"P| 2 = \T + + T~\ 2 , and |T cex | 2 = 2|T"| 2 , (3.2) 



in a non-trivial way such that they cancel in the isospin limit P 10 l In [6], the infrared divergences are 
removed at the amplitude level just replacing m 7 by the detector resolution E mSuX . As we have seen 
in Sects. 12.9.21 and 12.10. 2l this admittedly catches the leading effect, but neglects finite terms, which 
e.g. affect the logarithm such that effectively m 7 is replaced by 2P max . Hence, the full outcome of the 
soft-photon approximation should somehow be accounted for. For that reason, we define the partial 
waves in (I2.13P not for the pure scattering process T^n, but include bremsstrahlung via 

] o e 2 

T W N ~~ > TttAT+ttA^ = T^N + T wN ^g^ 2 ^'(' S ' ^maxji (3-3) 

where T^® is determined by (|2.10p and (|2.25p and C(s, t, E max ) is defined in (|2.77|) . The corresponding 
triangle relation is infrared finite, coincides with [7] at threshold and involves amplitudes satisfying 

|AW^v 7 | 2 = g Yl \T«N+wNy\ 2 = \M nN \ 2 + \M^\ 2 + 0(p 5 ), (3.4) 

spins 

such that the required decoherent summation of the pure process and bremsstrahlung is maintained 
up to higher orders. 

The problem we encounter for the P-wave projections is due to the threshold behavior of the 
partial wave amplitudes Tf. for elastic scattering, their real part can be generically written as 

ReI] = |p| 2Z [a x +k\p\ 2 + •••), (3.5) 

but once inelasticities arise, the prefactor has to be modified according to |p| 2 ' — > |p|'|p'|'. While the 
5"-wave is obviously not affected, the threshold behavior of the P-waves differs between the elastic 
channels and the charge exchange: whereas the former vanish ~ A (s,m 2 ,M 2 ), the latter behaves as 
A 1/2 (s, m 2 , M 2 ) A 1/2 (s, ml, M 2 ). Hence, R will be totally dominated by fi± (s) for s -> (m p + NU) 2 . 
Even more, if we take the denominator in the isospin limit, R will diverge as a square root at threshold. 

In order to circumvent this problem, we remove the known kinematical threshold behavior, i.e. 
instead of fi±(s) we consider the triangle relation for 

^ and SW, (3.6) 

P P P 



#10 As discussed in Sect. 12.101 we can even replace the amplitudes by their real part up to higher orders. 
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respectively. In this way, the quantity of interest for the P-wave projections becomes 

Rl±(s) ^y- , '^-^ l , ^\ ( 3,, 



which is not only regular at threshold !^ 11 ! but reduces just to a triangle relation for the P-wave 
scattering lengths, completely analogous to the .S'-wave projection. If we incorporate bremsstrahlung 
according to (|3.3p . (|3.ip (for I = 0) and (|3.7p should therefore yield a consistent framework to study 
isospin violation above threshold. 

3.2 Numerical results 
3.2.1 Low-energy constants 

In this section, we specify the required low-energy constants not yet collected in [7]. In the strong 
sector, these are C4, d\(/j,) + d v 2 (fJ.), d v 3 (fi), and d\ s (/x). The &x are often given in terms of the scale- 
independent quantities di, which are related to d\{n) by 

F^( M ) = F^ i + I ^log^ J for i€ {1,2,3}, i^f 18 ( M ) = F^d 18 + (3.8) 

C4 will be taken from [28], where various previous analyses are briefly reviewed yielding 

c 4 = 3.5+^GeV^ 1 . (3.9) 
dis is determined by the Goldberger-Treiman discrepancy 

A GT = 1- mp9A = - 2M ^\ (3.10) 

Ftt9ttNN 9 a 

There are many values for Aqt available in the literature [13,29-32], ranging from (1.4 ±0.9)% up to 
a few percent. However, small values of Aqt seem to be favored: first, the ex perim ental value of g& 



has increased over the years, whereas i*V decreased due to radiative corrections !^ 12 ! Both effects tend 
to reduce Aqt- Second, it is shown in [29] that only rather small values of Aqt ar e consistent with 
the standard picture of chiral symmetry breaking. We conclude that values of Aqt much larger than 
about 2 % seem very unlikely. Consequently, we will use [31] 

A GT = (1.4 ±0.9)% (3.11) 

in the numerical work, which corresponds to (// = 1 GeV) 

F%d\ 8 = (-2.1 ± 2.6) • 10~ 3 . (3.12) 

The remaining d\ are fixed based on ttN threshold parameters (for details of this method we refer 
to [33]): as far as LECs are concerned, the isovector S- wave scattering length a - , the isovector P- 
wave scattering lengths a^ ± , and the isovector S'-wave effective range 6q + , depend at 0(p 3 ) only on 
C4, d\ + d2, c?3, and d§. Inverting this linear system yields e.g. 

d +d = - 3 + 21 9a + 2 9a 9a (3 + 2gj) M n _ gj (4m 2 + 4m p M w + 3M 2 ) 

2 1152vr 2 F 2 1927rF 2 (m p + ik^) 32m 2 M 2 

2m p + M n ttF 2 _ _ _ . 

+ ™ vf- ( x — m p o 1 _+2m p a 1+ +3M,a 1+ . (3.13) 



* 11 In the actual evaluation, isospin violation in R will be restricted to first order in S, such that the correction factor 
in the denominator disappears. 

#12 Recently, F„ has further decreased from 92.4 MeV to 92.2 MeV, which is not included in the quoted literature. 
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For a and a 1± we use the same values as in [7]. Together with b 0+ = (8 ± 8) • 10 M w 3 this leads to 

F 2 {d\ + d\) = (38 ± 8) • 10~ 3 , F 2 d\ = (-38 ± 10) • HT 3 , F 2 d\ = (-2.4 ± 3.0) • 10~ 3 , 

c 4 = (3.4±0.3)GeV^ 1 . (3.14) 

Note that (|3.9h is reproduced quite well, while this value for d\ is at least compatible with [7] within 
the errors. For consistency reasons the new value for d\ is used in order to treat all d% in question in 
the same manner. Note that the errors of the di obtained in this way are strongly correlated. 

In general, the remaining electromagnetic LECs are treated as in [7]. However, there is one linear 
combination of LECs, for which this method fails, since the /3-function vanishes: 

e 2 F 2 {2g\ + g\) - 4B(m d - m u )(2d\ 7 - d\ 8 - 2d\ 9 ). (3.15) 

As except for d\ 8 none of these LECs is known, the central value of the whole combination is set to 
zero in the numerical analysis. One possible way to estimate the uncertainty might be based on the 
observation that this combination of LECs is nothing but g p — g Q (cf. (|2.30p ). In view of the finding 
for Agtj we regard 0.5% as a reasonable upper limit for the contribution of (|3.15p to isospin breaking 
in <7a- The same effect would be generated by a f3- function of the magnitude 

0.5 • 10~ 2 5A ^ ~ 11- (3.16) 

On the other hand, the vanishing of the /3-function means that this combination is particularly stable 
against variation of the renormalization scale, which hints at a value of the total counterterm not 
being exceptionally large. This is the reason why in [34] in such cases the /3-function is set to 1, which 
just reproduces the naive order-of-magnitude estimate 1/167T 2 . Since the /3-function of both 2g| + g\ 
and 2d\ 7 — d\ 8 — 2d\g vanishes, one could also argue that these two combinations should be treated 
separately. Counting e 2 F 2 ~ B(m& — m u ) and adding the errors linearly, this would correspond to an 
"effective" /5-function of ~ 5. All in all, we conclude that it is very difficult to obtain a reliable error 
for (|3.15p . What we will eventually use is 

e 2 Fl(2g\ + g\) - 4B(m d - m u ){2d\ 7 - d\ 8 - 2<f 19 )| ~ ^ ~ 3 • 10" 3 , (3.17) 
which should give a reasonable estimate, though not the most conservative possible. Note that this 



value is somewhat smaller than the estimate given in [35] in the context of the two-nucleon svstem P 13 



3.2.2 Threshold divergences 

In the numerical analysis, there are several obstacles to be overcome, which are related to the behavior 
of the amplitude at threshold. First of all, the basis functions of the tensor decomposition associated 
with Ig(s) (cf. diagram (sq)), Ig (s) (03), Vu(s) (t>i), and Au(s) (ag), and the box-grap h top ologies 



i] 3 (s,t) (^2), /7 3 (s,t) (07), Ai2(s,t) (as), and A2i(s,t) (a 8 ), involve prefactors of the tvp a^ 14 



1,1 , 

an d "T7 n 9 . , (3.18) 



A(s,m 2 ,M2) ' \{s,m 2 ,M%) + st' 



#13 The low-energy constants used in [35] can be most easily compared to the present framework by matching to the 
charge-symmetry breaking corrections to the pion-nucleon coupling constant given in [36]. 

#14 For the sake of simplicity, the following discussion refers to the isospin limit. However, we have checked that the 
conclusions persist if 7^ is taken into account. 
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respectively, which suggest threshold divergences even worse than the Coulomb pole. In fact, these 
singularities are canceled by the numerators, such that the loop functions are actually finite at thresh- 
old once the Coulomb pole is subtracted. To avoid any complications inherent to this (numerical) 
cancelation we start the analysis not directly at threshold sthr = (m p + M w ) 2 ~ 1.162 GeV 2 but at 
s = 1.17 GeV 2 . As the final results are smooth for s — > s t h r this implies no severe limitation. In ad- 
dition, the tensor decomposition of the box graphs causes difficulties also above threshold: it follows 
from (|2.5p that A(s, m 2 , M 2 ) + st = for z = —1 independent of the actual value of s, which induces 
numerical problems in the partial wave projection of the corresponding amplitudes. Fortunately, it 
turns out that the z-dependence of these loop functions is very weak, and, moreover, can be well ap- 
proximated by a straight line. For that reason, we replace the full z dependence by the interpolation 
between z\ = 0.8 and Z2 = 0.9. The precise choice of these points is of course rather arbitrary, but 
accounts for the fact that difficulties arise near z = — 1. 

Especially for the P-wave projections, another difficulty concerns the consistent treatment of the 
kinematics. Let us rewrite the relevant part of the numerator of (|3.7p as 

(E p + m p ){Af p (s) - Af p (s) - V2C A+ (s)A c r 

+ (v/i - mp) (Bf p (s) - Bf p (s) - V2C B+ (s)B^) } 
+ (E p - m p ){ - Al + ± l(s) + AfJ(s) + V2 C A ^s)A?^ 

+ (V~s + m p ) (^fW-^fW-VSCfl-W^)}, (3.19) 
where we have defined the correction factors 



C A ±(s) 



Cb±(s) 



(vg ± m n y - M 2 \( 8 ,m*,M*) 
\l (^i±m p ) 2 -M2 V A(s,m 2 ,M 2 )' 



(Vi ± m n ) 2 - M 2 ^ An \ j \{s,ml,M%) (g ^ 



\ (yfi± mp) 2 -Ml 2 (^5 T mp) J V A (s, m 2 , M 2 ) 



To avoid artificial threshold divergences, it is important to treat these correction factors in the same 
way as the corresponding part of the amplitude they are multiplied with. Therefore, we adopt the 
following conventions: for tree contributions (including the wave function renormalization) the full 
kinematics with An 7^ and A n 7^ are used. In particular, the nucleon masses in the axial parts of 
A LO and B LO in (|2.40p are restored according to (|2.3ip with Am = 0. Consequently, (|3.20p is applied 
as it stands. The sole exception is made for the counterterms involving e 2 ki or e 2 ^, for which the 
correction factors (|3.20p are set to 1, since the effects of the mass differences are anyway outweighed 
by the uncertainty in these LECs. In contrast, we have neglected the nucleon mass difference as far as 
loops and bremsstrahlung are concerned. Thus, for these contributions m n should be replaced by m p 
in (|3.20p . Note that to ensure consistency A„- is kept in the electromagnetic loops, although this is 
formally an C(e 4 ) effect. We conclude that the consistent treatment of the mass differences is crucial 
to prevent remnants of the different threshold behavior of the elastic channels and the charge exchange 
reaction as described in Sect. 13.11 from spoiling the threshold properties of R\±. 

3.2.3 Triangle relation 

The results for the triangle relation above threshold in the S- wave (R), the Pi- and P3-waves (Pol- 
and R\+), and the Q- and W-projections (Rg and Rn), are depicted in Figs. [T8H221 The central values 
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correspond to the LECs compiled in Sect. 13.2.11 a renormalization scale /j, = 1 GeV and a detector 
resolution E mSuX = lOMeV. The scale is varied between the mass of the p meson M p = 775.49 MeV 
and the A resonance ttia = 1.232 GeV, the detector resolution between 5 MeV and 20 MeV. The 
errors induced by the LECs are calculated via Gaussian error propagation, where the errors of the di 
are directly traced back to those of the corresponding threshold parameters. This ensures that the 
uncertainties in the different LECs are essentially uncorrelated. Finally, all three contributions are 
added in quadrature. We display the individual contributions to the triangle relation by tree graphs, 
strong loops, virtual photons, and bremsstrahlung, as well as the individual contributions to the final 
uncertainty. Due to the IR divergences, the separation of virtual photons and bremsstrahlung is scale 
dependent: by changing the scale in logm 7 //xrR, contributions can be shifted between virtual 
photons and bremsstrahlung. We choose //fr = m p , which emerges naturally in the calculation. 

Let us first consider the S'-wave. At threshold, the triangle relation is violated by about 2.5%, 
which is significantly larger than the 1.5% at threshold obtained in [7] (but within the uncertainty 
given there). The reason is that both tree and loop contributions have not been chirally expanded; 
in particular, the diagrams which at threshold only start at 0(p 4 ) are included. In this way, the 
difference between both values may be taken as an estimate of the potential impact of higher chiral 
orders. From this point of view, the error bands in Fig. [18] dominated by the uncertainty in the 
LECs would therefore not be modified substantially if the uncertainty due to higher orders were 
included. Remarkably, the residual dependence on the renormalization scale as well as on the detector 
resolution turns out to be negligible in comparison. The net effect of tree graphs and strong loops 
is rather small and practically independent of i/s. In contrast, virtual photons yield the dominant 
contribution, which grows significantly with increasing energy. For higher energies, bremsstrahlung 
becomes quite important: it counteracts the rise induced by virtual photons. The size of the finite 
pieces of bremsstrahlung apart from the leading terms oc log mp/4E^ ax amounts to about 20% of the 
total bremsstrahlung contribution. 

For the Pi -wave, we observe the same problem as in [5]: the denominator itself becomes very small, 
such that especially at threshold Ri- is very sensitive to the precise values of the LECs. We conclude 
that for this projection no meaningful statement can be made. In the other P-waves isospin breaking 
is very small; the triangle relation is violated by 1 % at most. The observation that the variation of 
Pmax yields a shift comparable to the uncertainty due to the LECs could be taken as a hint that the 
total error might be underestimated. Higher chiral orders and the A resonance might lead to larger 
shifts than indicated by the error bands in Figs. [T9l 1201 and!22i In particular the A resonance, which 
is only parametrized in the LECs of the theory and not included as explicit degree of freedom in the 
present work, could become important in the P-wave. Similarly to the S-wave we find that virtual 
photons and bremsstrahlung tend to evolve in the opposite direction for increasing energy, but at 
least for R-^ and Pi + these contributions no longer outweigh tree graphs and strong loops. Since the 
P-wave projections themselves vanish at threshold, bremsstrahlung can generate finite contributions 
in the full energy range, which is indeed observed in the case of Rg and P11 + . 



3.3 Comparison to earlier work 

Isospin violation above threshold has been the subject of earlier studies in the framework of heavy- 
baryon ChPT [5,6]. The findings in [5] differ quite substantially from ou r res ults: in the S- wave, the 



triangle relation is found to be violated by about —2.5 % instead of +2.5 % P 15 I and also for the P-waves 



(especially for the W-projection) significantly larger negative values for R are obtained. However, this 



* 15 Note that the scattering lengths given in the same paper yield R ~ 1 %, cf. [3,7]. The discrepancy is related to the 
question which threshold energies are inserted and again illustrates nicely that the final result for R is quite sensitive to 
mass differences. 
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Figure 18: Triangle relation in the S-wave. Upper panel: central result (solid line) with uncertainties 
(dashed). Middle panel: individual contributions to the central result. Lower panel: individual 
contributions to the uncertainty. 
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Figure 19: Triangle relation in the (/-projection. For notations, see Fig. 
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Figure 21: Triangle relation in the Pi-wave. For notations, see Fig. [THJ 
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analysis is incomplete in the sense that virtual photons were neglected. Indeed, we observe that in 
all cases their contributions shift the curves upwards, which partly explains the discrepancy. Further 
differences concern the definition of the isospin limit, the values of the LECs (the Cj and d{ are taken 
from [37]), and the treatment of higher chiral orders. The philosophy in [5] is to expand the amplitude 
strictly at 0(p 3 ) and only to include some higher-order terms whenever necessary to ensure the correct 
threshold behavior. As we have seen in the previous section, such higher-order effects can indeed lead 
to differences of ~ 1 %. Moreover, the denominator is not reduced to its isospin limit, which would 
generate in our work a dependence on some new combinations of electromagnetic LECs, involving 
e.g. k\ and k<i- As for the Pi -wave, the problem of the small denominator manifests itself in [5] in a 
different way: the curve is regular at threshold but diverges around yfs = 1.115 GeV. Since we can 
reproduce a similar picture by varying the LECs in a broader range than described in Sect. I3.2TT1 this 
does not point to a fundamental difference, but rather to the difficulties in the Pi-wave. Finally, we 
remark that our results agree with [5] in the qualitative energy dependence of the triangle relation in 
so far as R is found to be a monotonically increasing function of y/s. 

In [6] virtual photons are taken into account, but this work eventually focuses on strong isospin 
breaking, such that in particular electromagnetic contributions to the mass differences are switched 
off. A consequence of this point of view is that the mass of the ir°n state is larger than of ir~p, since all 
pion masses may be identified with the mass of the neutral pion once the tiny strong contribution to A n 
is neglected. Therefore, the three amplitudes entering the triangle relation no longer possess the same 
threshold. Nevertheless, we observe a similar discrepancy as in the case of [5]: the values for the triangle 
relation in the 5"-wave (about —0.7%), and in the Q- and especially W-projection, lie significantly lower 
than the present results. However, virtual photons (and higher-order contributions) again tend to shift 
the curves into the right direction. Although no reliable error bands are provided l^ 16 we therefore 



consider the origin of the deviations from former studies in ChPT to be essentially understood. 

It is very difficult to compare our results to phenomenological models [8, 9] , since in contrast to 
the analyses in ChPT an unambiguous separation of the different sources of isospin breaking cannot 
be achieved. In [8], a coupled-channel potential model [38] is used to remove Coulomb interactions 
and hadronic mass differences from the data. The triangle relation is found to be violated by about 
— 7% in the S'-wave at yfs = 1.10 GeV. This value is practically independent of -y/s (cf. Fig. 1 in [8]), 
such that a similar number could be deduced for the threshold extrapolation of the amplitude. In the 
"H-projection no isospin violation is found, whereas Rg displays an effect comparable to the 5-wave 



#17 



The investigation in [9] is based on a tree-level meson-exchange model [39] ; electromagnetic effects 
and hadronic mass differences are separated relying on the NORDITA method [40-42]. In the S- wave, 
R is found to be nearly independent of in the energy range ^/s = (1.10 — 1.14) GeV. However, the 
final value R = (—6.4 ± 1.4) % differs from the triangle relation for the scattering lengths obtained 
by extrapolating the amplitudes to threshold, which is violated by about —3.5%. The central values 
for Rg and R-h can be extracted from Fig. 5 in [9] and amount to (3 — 4) % and —(5 — 6) % in the 
considered energy range, respectively. However, the quoted uncertainties are quite large, such that at 
least in the W-projection the results are compatible with isospin symmetry. 

A similar analysis is performed in [10]. The corrections due to hadronic mass differences are 
calculated within the X-matrix formalism and found to be in reasonable agreement with the NORDITA 
results, while the electromagnetic effects are taken from NORDITA directly. For the S'-wave, isospin 
violation is found to be smaller than 1 % in magnitude. The discrepancy to [9] is traced back to 
the fact that a coupled-channel approach was chosen, arguing that elastic ir~p scattering and the 



* 16 In [5] only shifts due to the ignorance of /i and are considered, whereas the error bands displayed in [6] refer to 
the outcome of the MINUIT minimization routine of the CERN library, which, as pointed out by the authors, by no 
means reflect the total uncertainty. 

* 17 About —7% at yfs — 1.10 GeV. However, the situation is less clear at higher energies, cf. Fig. 2 in [8]. 
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charge exchange are coupled channels even if isospin is not conserved, whereas in [9] each reaction was 
analyzed separately. Isospin violation in the P- waves is not considered. 

We conclude that the phenomenological models do not provide a consistent picture. Moreover, the 
systematic uncertainties are very difficult to assess. As already remarked in [5], in principle a mean- 
ingful comparison between the model calculations and ChPT would require a detailed confrontation 
of the NORDITA method with our treatment of isospin- violating effects, such that the outcome of 
both methods cannot be naively compared. However, we have seen in the course of the numerical 
evaluation of the triangle relation that the final result (and especially the threshold behavior) is quite 
sensitive to the correct incorporation of mass differences and can in principle generate large effects. 
Certainly, these fine effects are not sufficiently well under control in the model calculations because 
they do not provide a systematic and consistent framework to analyze isospin breaking. 
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A Effective Lagrangians 

Our investigation of isospin-breaking effects relies on the effective Lagrangian for nucleons, pions, and 
virtual photons, as constructed in [20] 



r„ - r(p 2 ) + r^ 2 ) + Hp 4 ) + £( e V) , £ (p) , Ap 2 ) , A e ^ + + + r 

■'-eft — 7v ' 7r ' *-"k ' it ' '-N ' N ' N ' N ' N ' 



(A.l) 



Apart from the leading piece 



r(p 2 ) + r( e 

7T ' 7T 




'7 



^(d^uU^U + X ] U + U^ X ) + ZF 4 (QUQtf) - \f, v f^ - \{d^)\ 




(A.2) 



we actually need the following terms: 



4 eV) = F 2 ^(d fl U^d IJ ,U}(k 1 (Q 2 } + k 2 (QUQU^)) 



+ h((d^QU){d^QU) + (cPUQU^idnUQUi)) + k 4 {d"U^ QU){d^UQU^, 
4 p2) = *{ Cl < x+ > - £^(u,u u ) (D»D» + h.c.) + |<«X> + \*<^KM + C5X+ 




4 e2) = FH[h(Q\ - Q 2 _) + f 2 (Q+)Q + + h{Ql + Q- )}*, 

= *{ - ^(K- [D^]]D V + h.c.) - ^-(K, [D», u v \W + h.c.) 
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+ (t*V [D u ,u x ]](D»D»D x + sym) + h.c.) + d b — ([*_, u^D" - h.c.) 

+ d e ^([D^ F+]D» - h.c.) + <h±{[ir, (F+)]D» - h.c.) 

+ du (a^ {u u [D\ ,u^])D x - h.c.)+d 15 ^(a^( Ufl [D u ,u x })D x - h.c.) 

+ ^V75<X+K + ^r7 M 75(x+*V> + d ls l -r^[D^ x-] + di^^[D^ (x-)]}*, 
4 e2p) = F 2 *{^(Ql - Ql)^u, + |(Q + > V 75 n M + f 7 At 75(Q + )(Q+n M ) 

+ y7 M 75Q+(Q+^) + 96^({Q+)(Q-U^ - h.c.) + g 7 ^-{Q_(Q +Ull )D» - h.c.) 

+ gs^Q+iQ-ujD* " h.c.)}*, (A.3) 

where (A) denotes the trace of a matrix A, A = A- (A) /2 its traceless part, *(C+h.c.)* = *C*+h.c. 
for an operator O and 

d/JJ = d^U-iA^Q, U], x = 2.Bdiag(m u ,m d ), U = u 2 , 

= 8^A V - a v A„, Q = I diag(2, -1), J%, = T% = QF^ U , 

Q = ediag(l, 0), F% = (u)Qu ± uQu f ) , 

D li = d li + r M , r M = i («t(^ - iQA^)u + u{d^ - iQA„)u^ , 

X± = u^xu^ ± ux^u, = i(u\dfi - iQA^)u - u(d^ - iQA^u^, 

Q± = \{uQv) ± v)Qu), [Dp, u v ] = + [r M , u v \. (A.4) 

VP = (p, n) r contains the nucleon fields and the pion matrix U in the a representation reads 



TV 2 .T-TV ( 7T° V2 7T+\ 

W = \/l--^ + *— . r-7r=( ^ Q ). (A.5) 



F 2 F ' ly^Tr- 



-7T 



The renormalization of the LECs proceeds along the following pattern (the /3-functions are given 
in [20]) 

Pi Vi 
li = -/iX + l T i(fj,), h = Gi\ + kl(fi), (U = ^X + d-(fj,), gi = ^X + gKfj,). (A.6) 

Note that the choice of operator basis in corresponds to [17]. As discussed in detail in [20], 

changing this Lagrangian to the version used in [43] results in a redefinition of several LECs and their 
/3-functions. In particular, some dj are replaced by di according to 

d\(») = di(») + l jt0, d\ 8 (») = <rM-^L-ii(t,), 4(/i) = d^) + i |^). (A.7) 

Counting e ~ 0(p), the full calculation at order 0(p 3 ) involves terms of order 0{e A p~ v ), which are 
numerically tiny and were therefore dropped in [20] . Although we neglect these effects in the numerical 
analysis too, which reduces the number of LECs considerably, the amplitude shall be given in its full 
generality for the following reason: the cancelation of UV and IR divergences serves as a powerful 
check on the calculation, and the consistent treatment of C(e 4 ) terms ensures that this cancelation 
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works for all contributions at order 0(p 3 ). However, for this purpose, there are further contributions 
to the effective Lagrangian to be taken into account: first, a Lagrangian of the type C\ should be 
added to (lA.lj) , and second we need some terms which do not contain pion fields ( "contact terms" ) , but 
are indispensable in the context of the renormalization of the one photon exchange. Both pieces can 
be taken from [44], since there the same basis for Cf is used, such that the additional Lagrangians 
are consistent with [20]. Eventually, it turns out that we only need the counterterms encoded in 

^contact = -h^r* + Tl v T\ w ) + h,(Q 2 )F^F^. (A.8) 

The first part is already given in [17], whereas the second one is specific to the inclusion of photons, 
since the electromagnetic field strength tensor appears separately. As in [44] no counterterm is specified 
to absorb the pertinent UV divergence as occurring in the renormalization of the generating functional, 
we extend the operator list of terms without pion fields as set up in [17] by a forth term cx /14. The 
corresponding /3-functions read 

h i = 8 i \ + h\(n), S 2 = ^, ^4 = -^. (A.9) 



B Loop functions 
B.l Definitions 

In this chapter, we will give explicit expressions for the loop functions originating from the integration 
of scalar propagators and for the basis functions of the tensor decomposition. Working in d < 4 
dimensions, the UV divergences are captured in the pole term 

^^{th-^-'^v}- (B1) 

where \x is the renormalization scale and 7e = — r'(l) = 0.5772156649 . . . denotes the Euler-Mascheroni 
constant. To ease notation, also the variants 

C , = 2m^ 4 |^ I -i(log47r-7E + l)} = 327r 2 A + log^, (B.2) 
Q = 2M?~ 4 |— - - (log47r - 7E + 1) j = 32vr 2 A + log -±, i G {vr, tt }, 

are used. In contrast, IR divergences are regularized by a finite photon mass m 7 and expressed by 

2 22 

L = log^I, L i = log^ = L + log^, iG{vr,7r }. (B.3) 
m z M- 

Our results may be compared to calculations taming both UV and IR divergences in dimensional 
regularization by the replacement 



Cir = 2m 



d-4 



1 



IR ^ 



I (log 4tt - 7E + 1) }^ -L - 1. (B.4) 



d-4 

The propagators are implicitly understood to carry a small imaginary part according to 

m 2 — p 2 — > m 2 — p 2 — ie, e — > + . (B.5) 
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We will make use of the so-called dilogarithm or Spence's function 

log(l-t) 



Li 2 (x) 



dt- 



t 



(B.6) 



in terms of which some integrals involving virtual photons can be done semi-analytically. Furthermore, 
the Feynman parameter integrals 



f(t,m) 



dx 



Am 



1 - x(l 



h(t,m,Mi) = hf(t) 
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log ( 1 — x(l — . 



arsinh 



2m 



Mi 



l-x(l 



(B.7) 
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2m 2 
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arctan ■ 



-A(s,m 2 ,M?) 



m 2 +Al 2 - 



+ 7TI 



if A (s,m 2 ,M?) > 0, 
otherwise, 



if A (s,m 2 ,M?) >0, 




otherwise, 



appear in several loop functions # 18 l Some of the loop integrals involving virtual photons can be 
obtained from integrals which are already present in the isospin limit by sending the appropriate 
meson mass to zero. Therefore, i and j run over {tt, vr°, 7} in Sect. [R2] (with M 7 -> m 7 ). 

Whenever possible, i.e. if no infrared divergences prevent the dispersion integral from converging 
in the limit of vanishing photon mass, we give the loop functions both in terms of Feynman parameter 
integrals and based on dispersion relations. For this purpose, we also need the imaginary parts of the 
t-channel diagrams, which only emerge for positive t, while the formulae for the real part in Feynman 
parameter representation are only valid in the given form for physical values t < 0. The dispersion 
integrals in some cases require one subtraction, which is specified by the subtraction point Sq. In 
general, the expressions for -u-channel diagrams follow by substituting s «-> u and S <-> A, such that 
they are given explicitly only if modifications are necessary. 

As far as the charge exchange reaction is concerned, A„- can appear both via pion propagators and 
via the external kinematics. To obtain representations valid for all channels in question, we will make 
use of 

[ for n~p 

otherwise, 



7r°n, 



M r , 



(B. 



#18 The imaginary parts are not shown, since they are treated separately in the following sections. 
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As remarked in 13.2.21 we will also take care of in scalar loop functions involving photons, though 
this is an C(e 4 ) effect. 

B.2 Scalar loop functions 

We adopt the following notation: the infrared (regular) part of any integral over m meson and n baryon 
propagators is in general denoted by I mn {R mn ), the finite pieces as I mn = I m n\x^o (Rmn = Rmn\\-*o)- 

B.2.1 Meson integrals 

The result for the tadpole reads 

a -p - 1 t ddk 1 M ? r moi 
Ai = /l0 " 7 Jj (2^ Mi - h* - m^ Cl - (B - 9) 

For the integral over two meson propagators we find (Jij(t) = I^oi^)) 
1 f d d k 1 



Jij(t) 



i ( 27T ) d (M? - k 2 ) f Mf - (k - A) 2 ) " * J (t'-t)(t'-t o y 
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J ( + \± t - *° I dt'ImJ i3 (t>) 



1 } MfAx) — x(l — x)t 
h J (t) = - J&xlog ^ V m2 K - + M? j (x) = xMf + (l-x)Mj, 

o * 

M 2 -M 2 Mi W(M?,M?,t) Mi + Mf-t-\V*(Ml,Mf,t 
= — 2 H log — H — log 



2MiMj 



ku{t) = J dxlog - j^x(l - x)j = -2 + 2 J 4M [ t t arsinh 



2Mi 



eU-iMi + Mj) 2 ) 9(t-4M?) I AM 2 



hnMt) = -± — -A 1 / 2 (Mf,Mj,t) , ImJ M (t) = v — ' > yj 1 

J 77 ( t) = _ 2A + _i_| 1 _l g^|, i m J 77 ( t) = M. (B.10) 



B.2. 2 1 meson, 1 nucleon 
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B.2.3 2 mesons, 1 nucleon 

is symmetric in z and j. The regular part will be denoted by R l 2 J i(t)- 

21 * Ji ( 27T ) d (M 2 - k 2 ) [M] - (k - A) 2 ) (m 2 - (p - A;) 2 ) ' 
= ^2 + ^ + 2) /(t,m) + (tf(t) + , 

l2l{t) = 32^ (Q + 1} m) + 32^ ^ ( *' m ' Mi) + h2 ^ ' 



= + log ^ + l) /(*, m) + ^ Wit) + fcjr (t)) , 

1 i g ^-W-^ 1 -^ 
fc? (t) = i | dx^^ + (i - j), M 2 (x) = xM 2 + (1 - x)M 2 







2m 2 J (t-y/-t (4m 2 - t) \ (t + y/-t {Am 2 - t) 
^ (t) = V-t(4m 2 -t) l Ll2 1 2^ J " Ll2 I 2m~ 2 



(2m 2 - M 2 (x)) arccos -A^Wi-*)* 

h 2 {t) = 2 / dx =. 

(1 - x(l - x)4r) Jm 2 j(x) (4m 2 - M 2 -(x)) - 4m 2 x(l - x)t 

1 1 

h«(t) = 2 J dx- 



m M A a r ccns -M^sq-s)* 
Ml ~ 2^ J arccos 



!y(m 2 -i(l-i)()(M, 2 -a:(l-3;)() 



(i-xa-x)^) ^i-^-x(i-x) 



t 

i\i 2 



-^-x(l-x)t 

m arccos — - 



^rym = 2 f dx y m 2 -x(i~x)t 

2 J (l-x(l-x)^) V-x(l-x)t' 
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Sir y/ (Am 2 -t)t ^ ( t _ M 2 - M 2 ) 



V~t(t-Ml 


'-Ml 




hi- 


4m 2 )\^M 2 ,Mf,t^ 


2 \ 


/ m 2 X | 


(M?,M 2 ,t) 


i + M 2 M 2 t 



H = log 

87^ (i-4m 2 ) 



; . (t - 4M 2 4m 2 - t V (4m 2 - 1 (t - 4M 2 

Im/^ (t) = -i V ^ arctan - 

21 W 87r v / (4m 2 -t)t t-2M 2 



9 (t - 4m 2 ) t - 2M 2 + J(t - 4m 2 ) (/ - 4M 2 ) 

log 



8vr^t (t-4m 2 ) 2^/m 2 (t - 4M 2 ) + Mf 



8vr v /(4m 2 - t) i V i 8vr^t (t - 4m 2 ) 2m 

oo 

1 f diW^ y tj C + l J#(t) 

1 2l{t)~ n f t ,_ t «2lW, «2lW- 32vr 2 m 2A*'^ 327r 2 m 2: 

1 (2m 2 - M 2 ^)) arccos 2 ™ 2 - M y*> 



h$(t) = 2 J dx 



(1 - x(l - x)4r) A /M 2 -(x) (4m 2 - M 2 .(x)) - 4m 2 x(l - x)t 



1 / m Mi\ 2m 2 -Mf 

(l-x(l-x)^)./l-^-x(l-x)-^ 



: m arccos — 



hVit) = 2 / dx - ,ff7^ ^— ■ (B.12) 

B.2.4 1 meson, 2 nucleons 

The scalar loop integral over 1 meson and 2 nucleon propagators in its general form reads 

ri = i [ ^ 1 fB 13l 

12 i h [2n) d (M 2 - k 2 ) (m 2 - { Pl - k) 2 ) (m 2 - [p 2 - k) 2 ) ' K ' ' 

However, we will only consider two special cases relevant for our calculation, in which either both p\ 
and p2 are on-shell or at least one of them. Inspired by [13], we write Ii 2 = I a for p\ = p\ = m 2 , 
(pi - P2) 2 = t and I12 = Ib for e.g. p\ = s, p\ = m 2 , (pi - p 2 ) 2 = M 2 . In the first case, we find 

ri (f) = I f^L 1 

AK ' i h {2ix) d (M 2 - k 2 ) (m 2 - (p - A:) 2 ) (m 2 - (p> - k) 2 ) ' 



47 



* arccos 



, 2^Jm?-x(l-x)t 

gi{t) = / dx 



1 - x(l - x)-U J\-M^- X (l- X )J, 





t 2 \ + _ /l™2 i A/r2 



A 16ir v '(((-4m 2 ) M? AW Jr J f - f AW ' 

4m 2 

ImlVt) = \ ' log = L 

AK ' ltoy/t (t - Am 2 ) V m 2 

^k(t) = 32?r i m2 (C + 1) f(t, m) + f h(t, m, m) - ^(7 ) j . 



m 



arccos 1 



f ' 
/ dx 

I (l-x(l ■ - ' " 



= / dx (B.14) 



' m z / V 4m z v ' m z 



The definition 



d d k 

i J l (2vr) d (M 2 - A; 2 ) (m 2 - (pi - k) 2 ) (m 2 - (p 2 - A;) 2 ) ' 



^to = 7 / ^377777^ ,,,w,„, 77 7777779 7 7777 (B-15) 



for the second case needs some comment: the momenta are assumed to fulfil the on-shell conditions 
p\ = m 2 and p\ = s (or p\ = u for the u-channel). Apart from M, the masses of the external pions, 
which may enter through the kinematics, can contribute to isospin breaking. To account for this effect, 
we introduce a pion mass Mj that always corresponds to M n for the elastic channels, but has to be 
modified for the charge exchange reaction according to 



M 3 = », ., , . = S, (B.16) 

M^o lfpi = p' 



for s-channel diagrams and 



M i = < „, ., > V2 = A, (B.17) 
[ M T o if pi = p 

for the u-channel case. Mj enters ()B.15p by 

(p 2 -pi) 2 =M 2 . (B.18) 

The result for 1^ (s) reads 

#w = -32^5 (^) + x ) /"(-.^^) + 



9ij (8,X) 



fidxgV(s,x) if\{s,m 2 ,M 2 ) < 0, 

/^(s) ^^(s, + Jo b(S) d%9if(s, x) otherwise, 
2m 2 (M 2 + (1 - x) (s - m 2 )) arccos -Mf-(i-x)(s-m 2 ) 

V V /y 2M iy J(l-x)s+x(m 2 -(l-x)M 2 ) 

(1 - x)s + x (m 2 - (1 - x)M 2 ) } ^/-s 2 .(x) 
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„ , o , „^ M 2 +(l-x)(s-m 2 )- J s 2 .(x) 

2m 2 (Mf + n- x )(s- m 2 )) log V { ' ^ * 3 \ 

K 1 V " 2Mi^(l-x)s+x(m 2 -(l-x)Mj) 

{(1 - x)s + x (m 2 - (1 - x)M 2 ) } ^/4(x) 
s?-(x) = ((1 - x) (s - m 2 ) - M 2 ) 2 - 4M 2 (m 2 - x(l - x)M 2 ) , 
^ — — j (s - m 2 ) 2 - M 2 (s - m 2 + 2M 2 ) 

( a _ m 2)2 _ 4M 2 M 2 ^ V ; * V J J 



(?)( \ 

9ij (s,x) 



x l l{s) 



Im I l g (s) 



- 2M i] Jm 2 {s - m 2 ) 2 - M 2 M 2 (s + 3m 2 — M 2 - M 2 

1 9l (l-i)s+3)(m 2 -(l-3;)M, 2 ) 

/ m lo S m 2 (l-x) 2 

ff7j ( S ) = / dx 7 -y, 

J (1 - x)s + x (m 2 - (1 - x)M 2 J 

oo 

(m+M,) 2 

(m + Mi) 2 ) 
167rAV2 ( s ,m 2 ,M 2 ) 

| ( a _ m 2 + M 2) / s + m 2 _ M 2\ _ 2sM 2 + A l/2 m2) M 2) A l/2 ^ ^ M jA V 

X log— 

is \m 2 (s - m 2 ) 2 - M 2 M 2 (s + 3m 2 — M 2 — M 2 ) } 

ef s -(m + M^ 2 ) s ( s -m 2 - 2M 2 ) 
M j 16ttAV2 ( s , m 2 ,M 2 ) g m 2 s - {m 2 - M 2 ) ' 

61 f s _ m 2\ s + m 2 - M 2 + A 1 / 2 f s, m 2 , M 2 ) 
l m /^ (s) = ^f^J lpg { _K i± 

8ttAV2 ( s ,m 2 ,M 2 J 2m V^ 
^ (S) = 32^? (C+ l)/n( S ,m,M,) + — (M*) 



(M 2 \ 
1 - x(l - x)^- 1 

/ ux 7 

— i? A// 



/^(s) = h-yj{s) = / dx - v , 

(1 -x)s + x (m 2 - (1 - x)M 2 J 

_ (s)= . fite$\*,x) y ifA( S ,m 2 ,M 2 ) <0, 

[ dx ^ 1) ( s ' x ) + /o Cs(S) dx ^ij 2) otherwise, 

o/„,o , x / ow m 2 (l+x) + (l-x)(s-2xM 2 )-M 2 

2m 2 (M 2 + (1 - x) (s - m 2 )) arccos 2 — ^ lI l - 



9$ (.*>*) 



2j((l-x)s+x(m 2 -(l~x)M 2 ))(m 2 -x(l~x)M 2 ) 



~(2)/ \ 

9ij {s,x) 



{(1 - x)s + x (m 2 - (1 - x)M 2 ) } y/-4(x) 



m 2 (l+^)+(l-a ; )( S -2^Mj ! )-M2 + A /4(^) 

2m 2 (M 2 + (1 - x) (s - m 2 ) ) log / ' 

V * y " 2j((l-x)s+x(m 2 -(l-x)M 2 ))(m 2 -x(l-x)M 2 ) 

{(1 - x)s + x (m 2 - (1 - x)M 2 ) } ^/s?.(x) 
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i 9l (l-x)s+x(m 2 -(l~x)Mf) 

m log Tr, TT71 

m 2 ~x(l-x)M 2 



M*) = / dx — ^ -y. (B.19) 

j[ (1 - x)s + x (m 2 - (1 - x)Mj) 

B.2.5 1 meson, 3 nucleons 

The so-called box graph leads to the scalar integral 

1 f d d k 1 

*L3(*. t) = J ^ ^ M 2 _ fc2 ) ( m 2 _ (p _ fc)2) ( m 2 _ (VJ _ fc )2) ( m 2 _ ( p / _ fc)2) ' (B - 20) 

At least the imaginary part can be given in a rather compact form 

9 (s - (m + M^ 2 ) 



Im/ 13 ( S ,t) 



2Ci(s) - iA (s,m 2 , M 2 ) + A ccx C 2 (s) + J-tA (s, m 2 , M 2 ) 
x log 



2v / Ci(s)|m.^McoxCi(s) 

6(s-(m + Mi) 

lmP 13 (s,t) 



e(s-{m + Mi) 2 ) J-tX (s,m 2 ,M 2 ) 

arsinh ■ 



- x=0 4^^^/4(1(5) -tX (s,m 2 ,M 2 ) VCi(s) 

= m 2 (s - m 2 f - M 2 M 2 (s + 3m 3 - M 2 - M 2 ) , 
( 2 (s) = M 2 (s + 3m 2 — M 2 — 2M 2 ) , 

A 2 



C 3 (s) = 4Ci(s) - tX {s,m 2 ,M 2 ) + 2A ccx C 2 (s) + M 2 (4m 2 - M 2 ) 

,2\ 



(7(5 — TTl 

hnl7Js,t) = i I - ^inh 1—1 (B.21) 

13 4vr (s - m 2 ) ^/=t (4m 2 - i) 2m 

After performing the integration corresponding to the extraction of the infrared part, we are left with 
an integral over two ordinary Feynman parameters. However, the representation 

1 1 

I \z{s,t) = j dx J dyj^f 13 (s,t,x,y), 


M 2 + (1 - y) (s - m 2 ) 



f[ 3 (s,t,x,y) = 



s\ 3 (s, t, x, y) {4M 2 4( S , t, x, y) - (M 2 + (1 - y) (s - m 2 )) 2 } 

,,^-9 -M 2 -(l-y)(s-m 2 ) 

AM 2 arccos — — >- 



2M iy J s\ 3 {s,t,x,y) 



,9-1 3/2' 

{4M 2 4 3 ( S , t, x, y) - (M 2 + (1 - y) (s - m 2 )) 2 } 
s\ 3 (s, t, x, y) = (1 - y)s + ym 2 - x{\ - x)y 2 t - y(l - y) (M 2 - (1 - x)A ccx ) , (B.22) 
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is numerically only viable below threshold. One possibility to obtain a representation valid above 
threshold is to employ dispersion relations 

oo 11 

4(M) = ~ J ds ' Im s> I [ 3( s S,,t) -R\ 3 (s,t), R\ 3 (s,t) = jdxj dy^fUs,t,x,y), 

(m+Mi) 2 

AMf arccos — - - '- 



fUs,t,x,y) = — — ^ 

{4M?s\ 3 (s,t,x, y) - ((1 -y)(s- m 2 ) + M 2 ) 2 } 

^ £-y)2( s - m 2 ) 2 + M?(l-y)(s-m 2 )-2M?s\ 3 (s,t,x,y) 
s\ 3 (s,t,x,y) {s\ 3 (s,t,x,y) - (1 - y) (s - m 2 )} 
1 

X 4M 2 s\ 3 (s, t, x, y) - ((1 -y){s- m 2 ) + M 2 f ' 



(B.23) 



where we have discarded the t-channel cut, which lies outside the low-energy region. The regular part 
R\ 3 (s, t) is numerically uncritical. Unfortunately, (|B.23P does not apply to Ij 3 (s, t), since Im Ii 3 (s, (X 
1/ (s — m 2 ), such that there is a non-integrable singularity at s' = m 2 . However, this divergence 
should be captured by a finite m 7 leading to a term proportional to L. Thus, the dispersion integral 
is certainly not suited to extract the infrared divergence. As described in Sect. I2T8TT] a representation 
appropriate for this purpose can be found by performing the y integration in (|B.22|) , Below threshold, 
the result is given by 



i 



^isM) = J dxg{ 3 (s,t,x), 



o 

,2 



i , . s 1 \s-m 2 2/1 (1- 2/2) A 1 + y A 2 , c ly /c^ 

£13 ( s > <> x ) = \ lo § / TT 1= arctan 

cic 2 [2/2-2/1 -2/2 (2/1 - 1J yC3 s - C12/0 

-42 , (1 - 2/o) 2 + c 3 -Bi - A22/2 , 1 - 2/2 , -Bi - -422/1 , 2/1 -1 

2 2/o + c 3 2/2 - 2/1 -2/2 2/2 - 2/1 2/1 

8M 2 ( 2/1+2/2-22/12/2 -Mi 

H — tt^ S — , arccos ■ 



(2/1 - y 2 ) 2 1 - 1) (1 - 2/2) 2^/(1 - 2/0) 2 + c 3 



3/2 

c 2 vyi — y2; ^ \ 

n . m 2 — s — M 2 1 

-2V=^axccos I 

+ ,„ 4M - ( 2d, log ^lll^L + ^+2/oA arctan . Cl ^ 



C2 /2 (2/1 - 2/2) 2 V 7 ^ I yi ~ 2/2 -2/2 (2/1-1) ^3" S - C12/0 

A 4 (l-2/o) 2 + c 3 £ 2 + 2/1^4 , 2/1-1 5 2 + 2/ 2 -4 4 1 - y 2 \ 
2 2/o + c 3 2/2 - 2/1 2/1 2/2 - 2/1 -2/2 J 

Ai(s, t, x) = ~^^{ (s ~ m 2 ) (ca + 2/0 " 2/12/2) + (a - m 2 + M 2 ) (2/1 + 2/2 - 2^) }, 

A 2 (s, i, x) = — I (s - m 2 ) (c 3 + 2/0) (2/1 + 2/2) - 2 (s - m 2 ) 2/02/12/2 



(s - m 2 + M 2 ) (c 3 + 2/0 - 2/12/2) }, 



51 



A 3 (s,t,x) = ( c 3 + 2/0 2 ) 2 (2/1 + 2/2) + (03 + 2/0) 2/12/2 (yi + 2/2 - %o y 

+ d 2 ( (03 + 2/o) ( 2 2/o (2/1 + 2/2) - y\ - yi) + 42/o2/i2/2 (yi + 2/2 - 2y ) - 2y?y|) }, 
A 4 (s, t, x) = — jdi ( (c 3 + y 2 ) (yi 2 + 2/2 2 ) + 2yi2/ 2 (2/1I/2 - 2/0 (2/1 + 2/2)) ) 

+ cfe ( (03 + 2/o) (2/1 + 2/2) + 2/12/2 (2/1 + 2/2 - 4yo) J } , 
Bl(s,t.aO = ^{ (s-m 2 +Mf) yi y 2 ( yi + y 2 -2y ) + (s - m 2 ) ((c 3 + y 2 ) 2 

+ ( c 3 + 2/o) (2/1 + 2/12/2 + 2/1) - 2 (03 + 2/o) 2/0 (2/1 + 2/2) - 2yo2/i2/2 (2/1 + 2/2 - 2y ) J } ; 
B 2 (s,t,x) = -Jr{ d i{ ( c 3 + 2/o) (2/1 +2/2) + (c 3 +2/o) (2/1 + 2/2) 3 - 2 (c 3 + y 2 ) y (2/1 + 2/2) 2 

- ( c 3 + 2/o) (2/1 + 2/2) 2/12/2 + 22/12/2 (2/12/2 - 2/0 (2/1 + 2/2)) (2/1 + 2/2 - 2yo) ) 

+ d 2 yW2 (2 (c 3 + 2/0) + 2/1 + 2/2 - 2 2/o (2/1 + 2/2) ) } , 
iV(s, t, x) = (c 3 + (yo - 2/i) 2 ) (c3 + (2/0 - 2/2) 2 ) , 

di(s,t,x) = M 2 + (s - m 2 ) (1 - y ) , d 2 (s,t,x) = (s - m 2 ) c 3 - 2/0^1, 
d 3 (s,t,x) = (s - m 2 ) 2 - 4M 2 ci, 

c 4 ± ^/c 2 -c 2 A(s,m 2 ,M 2 ) s _ m 2 + M 2 _ ^ _ x ) Accx 



2/1/2(5, *, a?) = , 2/o(s, t, x) 



c 2 ' " uv ' ' ' 2(M2-x(l-x)t-(l-x)A cex )' 

ci(/,x) = M 2 - x(l - x)i - (1 - x)A cex , c 3 (s,i,x) = — - y 2 , 

ci 

c 2 (s, t, x) = (s - m 2 f - 4M 2 M 2 + 4x(l - x)M 2 i + 4M 2 (1 - x)A ccx , 

c 4 (s,t,x) = (s - m 2 ) (s-m 2 - M 2 ) - 2M 2 M 2 + 2M 2 (1 - x)A cex , (B.24) 

where, for convenience, we have omitted the arguments on the right hand side of the equations and 
suppressed the index i in the case of the auxiliary functions. Note that below threshold, the kinematical 
region where 



m 2 - mJ 4M| - (t ^ ccx)2 < s < m 2 + mJ 4M% - (t ^ cex)2 (B.25) 

has to be be paid attention to, since c 2 (s, t, x) can become negative and y\/ 2 (s, t, x) complex. However, 
the representation (1B.24|1 can still be used in this range if the analytic continuation is done properly, 
which can be achieved by evaluating the amplitude at s — > s + ie. As we are mainly interested 
in the continuation above threshold, we will not pursue this any further. While y 2 (s,t,x) < for 
A (s, m 2 , M 2 ) < 0, we have < y 2 (s, t,x) < 1 for A (s, m 2 ,M 2 ) > (and yi(s, t, x) > 1 in both cases). 
Taking this into account, the continuation of (IB.24jl above threshold reads 

1 

P 13 (s,t) = J dxg\ 3 (s,t,x), (B.26) 


a* (six)- 1 { s ~ m \ oz y^ 1 ~ Al + ^ 2 arctan 



cic 2 I 2/2 - 2/1 2/2 (2/1 - 1) y 7 ^ s - ciy 
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(1 - Uo) + c 3 Bi - A 2 y 2 , \-y 2 B x - A 2 y x yi~l 

2 2/o + c 3 2/2 - 2/1 2/2 2/2 - 2/1 2/1 

8M 2 J yi + y 2 -2yiy 2 -M< 

H — tt^ S — , arccos 



C 2 



;/2 (2/1 - V2? I V(yi-i) (1-2/2) " " 2 ^r^ (1 _ yo) 2 + c 3 

s - m 2 + Mf - Jx (s,m 2 ,M 2 ) 



+ y/mm log 



4M ^ J oj 2/1 + 2/2 , 2/1 (1 " 2/2) . A3 + 2/0-44 , ci^3 
+ TJo — o ^ 2«i log — — H — arctan ■ 



c 2 /2 (2/1 - 2/2) 2 1 2/1-2/2 2/2 (2/1 - 1) ^ a - ciy 

, At (l-2/o) 2 + c 3 , ^2 + 2/1-44, 2/1-1 #2 + 2/2^4, I-2/2] 

+ "V lo g 2~ + lo S lo § ( • 

2 2/o + c 3 2/2 - 2/1 2/1 2/2 - 2/1 2/2 J 

Based on ([R26]) . it is now straightforward to perform the limit m 7 — > 



8vr 2 (s - m 2 ) y^-t (4m 2 - t) 1 V«- m V J 16vr 2 



/fs + m 2 — M 2 + (1 — x)A cex ) arctan — — 2 — rr— t -\ — 
dx y ; \ u \ \, — s +m ^ H ,- x) ^ {B27) 
[s — m z ) [m A — x(l — x)t) \/s{ 3 

sj 3 = 4s (M 2 - x{l - x)t - (1 - x)A cex ) -( s -m 2 + M 2 -(l- x)A cex ) 2 . 
B.2.6 1 meson, 1 photon 

1 f d d k 1 = _c JL -i 

10 i J l (2vr) d (m 2 - k 2 ) (M 2 - (q - k) 2 ) 16vr 2 ' 1 ' ' 



B.2.7 2 mesons, 1 photon 



()= l fd^k 1 = 1 ] a^mMO 

20U » A (2vr) d (m 2 - A; 2 ) (M 2 - (q - k) 2 ) (M 2 - (q> - k) 2 ) vr J t'-t 

Ln f( h{t,M nj M n ) 
>J V'j 7r ) ~~r 



32n 2 M% J v ' w; 32vr 2 M 2 
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B.2.8 1 meson, 2 photons 



10 



( ) = 1 f d d k 1 , ^ _ t) 

U i J j (2vr) d (m 2 - A: 2 ) (m 2 - (jfe - A) 2 ) (M 2 - (g + A:) 2 ) {Q ~ q) 

1 r d d k 1 _ 1 fdt'lmP 10 (t') 

Jl (2ir) d (m 2 -(a- k) 2 ) (m 2 - (q 1 - k) 2 ) (Ml - k 2 ) ~ nj ' 



i J l (2vr) d (m 2 - (q - k) 2 ) (m 2 - (q' - k) 2 ) (M 2 - k 2 ) irj t'-t 
32^^ 1Og i^ + 16^M| + /ll0(t) ' 



1 

(*) = y r 



dx log a; 



o 



x(l — x) 



Mi 



M 2 j L . / 1 - y-t (4M 2 - g \ L ^ ( t + y-t (4M 2 - t) 



y/-t(AM 2 -t)\ \ 2M 2 I \ 2M 2 



h w (t) 



i 

M n f dx 



32tt7 (M 2 - x(l - x)t) y/-x(l - x)t l^yj-t (4M 2 - t) 



Im Plo(t) _ mO(^-t) arctan fiW7 1 + ^-4^) aKosh 



8vr^(4M 2 -t)t V t 8tt0 (* - 4M 2 ) ' ' 2M* 

Pio(i) is merely needed in the tensor decomposition of Pu(s,t). 

B.2.9 1 meson, 1 nucleon, 1 photon 

F M _ 1 f ddk I = ( / i\ 

11 {S> i J l (2ir) d (m 2 — k 2 ) (m 2 — (p— k) 2 ) (M 2 — (q + k) 2 ) W ~" 9 ' P ~* P j 

{/in(s) - L^n (s,m,M^-)} - i?n(s), 



32vr 2 m 2 



C + 1 1 

#11 (*) = ~ QO 2 2 /n (s,m,Mi) - fen (5), 



'11 



(-)=/ 



J. log( X+ (l- X )^|- X (l- a; )-^) 
<Lr 



x+(l-x)^f -x(l-x)-* 



, /x "i 2 I T /2s — st, \ /s+\ 1, , 2s -at, 

^1 s = TT727 <T17<d 2Ll2 o 11 " 2Ll2 — + 9 lo § o " 

A 1 ^ ( s ,m 2 ,M 2 ) [ \2s-s u J \s u J 2 2s - s n 

1 , s+ X(s,m 2 ,M 2 ) (2s- sti) sT-,) , n9 

- 2 ^ 2 ^ + 108 los (2.-.r!).a ) ,f - £ (m - M *> • 

2ir 2 m 2 m 2 j (2s-s^\ f s Ti\ 2tt 

hll{s) ~ W 2 (s,m 2 ,MZ) + \V 2 (s,m 2 ,M%)\ 2Ll2 V^^J + 2 Ui J " "3 

+ I log 2 2g ~ g n + I log 2 £11 + log A(g,m 2 ,M 2 ) (2s - a+) s^ \ 

2 2s - s^ 2 s n m 2 s (2s - sj^) s^ J 

if s > (m + M^) 2 , sf^s) = s - to 2 + M 2 ± A 1/2 (s, m 2 , M 2 ) , 
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i / m su(s,x)+x\ 1 / 2 (s,m 2 ,M%) 

\. sn(s,x log , ^ '= 

ku{s)= dx y-*- -ppr- if A (s,m ,M ) > 0, 

{ x \^( s ,m 2 ,M 2 )(l-x + x^-x(l-x)^) 

1 i \ sii(s.x) 

i siils, x) arccos , „ 

fcn(s) = dx — -/ " " M 2X ifA(s,m 2 ,M 2 <0, 

^ xV-A(s,m2,M2) (l - x + x^ - x(l - x)^j 

su(s,x) = 2m 2 + x(s -m 2 - Ml) , 

oo 

y 11 (s) = y 11 (s ) + i ? n(so) + ^ / -4^M^ -i. ll( s) ; 

7T J (S' - S)(S' - S ) 

(m+M n ) 2 

e( s -( m + M n ) 2 ) / A (s,m 2 ,M 2 ) \ / 
ImFn s = V ; ==L ( log ' — — LI. B.31 

For s > (m + M,,-) 2 , Fn(s) contributes to the Coulomb pole at threshold. The divergence is contained 
in hn(s) and solely due to the term oc 2ir 2 , while the remaining bracket is regular at threshold. Besides 
Vu(s), which emerges from vector-type diagrams, there is a second kind of loop integrals including 
one meson, nucleon, and photon propagator, arising from axial-type topologies 

1 f d d k 1 

ll(Sj ~ i J l (2vr) d (m 2 - k 2 ) (m 2 - (E - k) 2 ) (M 2 - (q - k) 2 ) ~ [q ^ q ' 



i 

dx- 



logx 

Ml 



x + (1 — x)-^- — x(l — X 

v ' m* v 



2 

m 



Li 2 



s - m 2 + M 2 - A 1/2 (s,m 2 ,M 2 ) 



AV2 ( s ,m 2 ,M 2 ) | y 2s 
' s — m 2 + M 2 + A 1 / 2 (s,m 2 ,M 2 ) 

~2s~ 



-Li 2 — -r^Y ifA (s,m 2 ,M 2 )>0, 



JM _ ^ g -L n -L _ ^ M 2 \ W ± <lZfj (^"' 2 )±( 2 *- ax (s.s) 



{(1 - x)s + (1 + x)m 2 - xM 2 } log ^ , 

dx- 



(x + (1 - X)^ - X(l - x)^f ) V*?lax(*,3) 

'. {(1 - x)s + (1 + x)m 2 - xM 2 } arccos 2 )-^-D^ 

f , LV ; V ; n} 2xM 7T J(l-x)s+xm 2 -x(l-x)M 2 
dx -; -r 

; , (* + (1 - ~ 1 " ) V^lax^) 

2 

xff(s) = k ^ if s ^ m 2 , 

11V ' s - m 2 + M 2 ± 2mM n < 

s 2 lax (s, x) = (1 - x) 2 (s - m 2 ) 2 + x 2 M^ - 2xM 2 ((1 - x)s + (3x - l)m 2 ) , 

oo 

a /x . , \ nax/ . s — sof ds'lmA n (s') „ x/ , 
A 11 (s) = An (so) + i?^(s ) + (g/ _ g)(s/ _; o) - ^(s), 
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9(s-m 2 ) s -m 2 + M 2 + A 1 / 2 (s,m 2 ,M 2 ) 
lmAll{s) = 8^m^W) l ° g ^M-fs ' 

^ (S) = -S^2 /ll(s ' m ' M7r) - 16^^ (S) ' (B - 32) 



{(l-x)s + (l + x)m 2 -xM 2 }log 



(l-x)s+{l+x)m 2 -xMl + yJ s 2 ^ ^(s,x) 



, , , 2mJ (l-x)s+xm 2 -x(l-x)M 2 

hff(s) = / dx v 



x 



o V m m 



i {(1 - g) a + (1 + x)m 2 - xM 2 } arccos 

f , iK ' y ' w > 2mJ(l-x)s+xm 2 -x(l-x)M 2 
OX -. -r 

(X + (1 - x)^ - x(l - X)^) V-«?lax(*^) 



B.2.10 2 mesons, 1 nucleon, 1 photon 

1 /* d d k 1 



i 7i (2vr) d (m 2 - fc 2 ) (m 2 - (S - k) 2 ) (M 2 - (q - k) 2 ) (M 2 - (q' - k) 2 ) 
f(t,M n ) \ ( m 2 V\ g 21 (t) 



A 21 (s,t) 



32vr 2 M 2 (s - m 2 ) \ & \s-m 2 J j 16vr 2 (s - m 

521 (t) = / dx ? 

/ (l - x(l - x)^) V4m 2 (M 2 - x(l - x)t) - M4 

j, i r 4m 2 arccos (i-y)«+(i+y)™ 2 -yMjj 
R 21 (s,t) = J dx J dy^) W3i(^ 



16vr 2 \ " ( s 2 i(s , t)X)y) )3/2 

(1 - y)s + (1 + y)m 2 - yM 2 } 
+ " i f 'J ) 7-—hr ' *' X ' 55 S ^ t, x, y)\ Accx=0 , 



S2i( s >*, z,y)s2i (s,i,z, y) 
s 2 .^, f , x, y) = 4y 2 m 2 (M 2 - x(l - x)t) - ((1 - y) (s - m 2 ) - yM 2 ) 2 , 

lmA 2 i(s,t) = i — - ' = arsinh . B.33) 

4vr(s-m 2 ) v/-t(4M 2 -t) 2M n 

We cannot write down the usual dispersive representation due to the non-integrable singularity at 
s = m 2 of the imaginary part, which, once regularized by a finite photon mass, generates the term 
proportional to L. 

B.2.11 1 meson, 2 nucleons, 1 photon 

Ai 2 (s,t) possesses both s- and it-channel cuts, is infrared divergent and contributes to the Coulomb 
pole at threshold. Note that we still consider s and t as independent variables and u to be fixed by 
the Mandelstam relation s + t + u = 2m 2 + 2M 2 — A cex . The s-channel integrals 

A r x 1 [ d d k 1 
A 12 (s,t) = - / — 
i J j (2tt 



1 f d a , 



) d (m 2 - k 2 ) {m 2 - (p - A;) 2 ) (m 2 - (S - fe) 2 ) (M 2 - (g' - fc) 2 ) 

1 (B.34) 



)d ( m 2 _ fc 2) ( m 2 _ ( p / _ fc )2) ( m 2 _ ( S _ fc)2) ( M 2 _ ( g _ ^2) 
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as well as 



d d k 



A 12 (u, t) . ^ ^ d ^ _ ^ ^ 2 _ ^ _ ^ ^ 2 _ ^ _ ^ ^ _ ^ + ^ (B.35) 

are only relevant for the elastic channels, while 

_ 1 /" d d /c 1 

12 («, *) - 7 ^ ( 27r )d ( m 2 _ fc 2) ( m 2 _ (p _ fc )2) ( m 2 _ ( A _ fc )2) ( M 2 _ ( g + fc )2) (j - >(>) 

does contribute to the charge exchange reaction. We will give explicit expressions for the s-channel 
case, where A cex is incorporated such that (|B.36P is reproduced correctly when deducing the u-channel 
via s <-> u. The result for the imaginary part reads 

Im^i 2 (s,t) = lmAl 2 (s,t) + Imi^s, t), 

9(s- m 2 ) e(u-(m + M,,-) 2 J m 2 + M 2 _ u 

Im AU(s, t) = v . ' arccos — -f 

87r(s-m 2 ) V-AKm^M^) 2mM,,. 

6> urn + M n ) 2 — u) m 2 + M 2 -u 

H , arcosh ■ 



87r(s-m 2 ) ^X(u,m 2 ,M 2 ) 2mM 7T 

0(u-(m + M n ) 2 ) r (m 2 -s) 2 A 2 (u,m 2 ,M 2 ) 
12{ , ) 16ir(m 2 -s)X 1 / 2 (u,m 2 ,M%)\ m 2 (A (u, m 2 , M 2 ) — M 2 A cex ) £(s, t) 

A CSL =o e(u-(m + M n ) 2 ) j (m 2 -s) 2 X(u,m 2 ,M 2 ) \ 

log 577 — 7; ^ ( 1 



16^ (m 2 -s)X 1 / 2 {u,m 2 ,M 2 )[ fa m 2 £{s,t) 
£(s, t) = 2M 2 [M* - A 2 CX ) - M 2 (M 2 - A ccx ) (3m 2 + u + 3A CCX ) + m 2 (u - m 2 ) 2 

Acc ^ =0 2Ml - (3m 2 + u)+m 2 (u- m 2 ) 2 . (B.37) 

Again, lmAl 2 (s,t) oc l/(s — m 2 ) prevents the direct application of dispersion relations. The real part 
is found to be 

#11 (it, m, M n ) J / m 2 \ 2 1 5i2(s, t) 

^ t} = 32vr 2 m 2 ( a - m 2 ) \ L + bg J } + 32vr 2 ( g -m 2 ) -^ 12(g ^ ) ' 

t> ( ,\ _ f i /, 1 f 2x 2 M 2 -gi 2 (x,y) 

12^, J y y i6 ^2j s 2 2(X)y)(m 2_ y(1 _ y)(M 2_ Accx)+x 2 M 2_ 5l2(X)y)) 







2sn{x, y) arccos 2(m'- 1/ (i- 1/ )(^-A cmt ))- 5ia (x, y ) 



+ 



2Vm 2 -s/(l-y)(A/|-A ccx )Vm 2 -^(l-y)(M|-A ccx )+x 2 M|-S 12 (x,j/) 



(s 2 2 (x,y)) 3/2 

s 2 2 (s, i, y) = 4x 2 M 2 (m 2 - y(l - y) {M 2 - A ccx )) - s\ 2 {x, y), 

s u (s, t, x, y) = (x- y) (s - m 2 - xM 2 ) + x(l - y) (t - M 2 + A cex ) + :r 2 M 2 . 
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For u < (m + M n ) 2 , gi2(s, t) is given by 



j*i2 («.*) dxg W ( S) t> x ) + /*i 2 ( ( *^ 6xg$ (s,t,x) + £+ (Si<) dx 5 f 2 } (s, t, x) otherwise, 



sit = ^{ M n (4m 2 + A ccx ) - M 4 + 2m 2 (m 2 - t - A ccx ) 

- i\W(M2 - A ccx ) (4m 2 - Ml + A ccx ) (4m 2 - M^)} Ac ^ =0 m 2 - t, 
± 1 , 4m7eCM) T 2 (2m 2 - M 2 ) A ccx T A 2 CX 

X 12V S ' r 7 — o ± 



2 



8 (m 2 u - (m 2 - M 2 ) 2 ) - 8A CCX (m 2 — M 2 ) - 2A 2 CX ' 



g$(s,t,x) 



2 ((1 - x) (2m 2 + A ccx ) + (2x - 1)M 2 ) arctan j^ ^^f^l^ 



jxM 2 + (1 — x)m 2 — x(l — x)u|y-uS 1 2 ) (s,t,x) 



|(l-a;)(2m 2 +A ccx )+M2(2a;-l)+^/«^ ) (s,t,a;)| 



(2), + , ((1 - x) (2m 2 + A ccx ) + (2x - 1)M 2 ) log W^+Ci-^^-^i-^n) 

g\ 2 '(s,t,x) = - 



u 



jxM 2 + (1 — x)m 2 — x(l — x)u^\J^^s^x) 

(s, t, x) = 4m 2 x(l - x) (u - m 2 - 2M 2 + A ccx ) - (4m 2 - M 2 ) M 2 (l - 2x) 2 

+ 2A CCX (2m 2 - Ml) (1 - x)(l - 2x) + (1 - x) 2 A 2 cx . (B.38) 



Above the u-channel threshold, we perform an integration by parts to improve the numerical properties 
of (712(5, i). In particular, the contribution to the Coulomb pole is separated into the boundary terms. 
It is solely given by the term proportional to tt 2 /\/ A (u, m 2 , M 2 ). Suppressing the dependence on s 
and t of the auxiliary functions on the right hand side of the equations, we arrive at 

(i) + 

T (2) 
x 12 



l 12 



dx 



(4 ) (x)4 1 2 )(x)-feg ) (x)/g ) (x)) - 



A 1 / 2 (u,m 2 ,M 2 ) 



+ ^ 1 2 ) (x)/ l g(x)-fcg ) (x)/ l g ) (x) 



„(2) 
'12 



,(1) 



x r2 M) - 2 ( Xl2 + 2^) ' x i2( s '*) - 2 (^i 2 + 2?) ' 



(a, t) = u + m 2 - Ml ± A 1 / 2 («, m 2 , M 2 ) , xf 2 < ^ < ^ < x+, 



2u 2u 



,„, . 2 ((! - x) (2m' + Aej + (2, - l)jg) log O-IC-^H^- W^V^M 
h 12 (s,t, x) — ^ 

V U 12 0*0 
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hf 2 \s,t,x) 



^12 ( S ' ^' X ) 



(1 - x) (2m 2 + A cex ) + (2x - 1)M 2 



2„,(2)^ (l-x)(2m 2 +A cex )+(2o;-l)M^+ A /M< 1 2 ) (x) 



dx / Q1, N \ 3 / 2 



« («12 fa)) 

(1 - x) (2m 2 + A ccx ) + (2s - 1)M 2 
u{xM 2 + (1 — x)m 2 — x(l — x)u} Ui2(x) 



((m 2 - Ml + (1 - 2x)u) V^SW + 4? (a;)) , 



(s, t, x) =- ^ (X) - 2m2u ^ ^ 



dx ( m , A 3 / 2 ' 



(u^O*))' 



2 



12 V ' ' ^ 2A 1 / 2 ( u , m 2 ,M 2 )' dx xM 2 + (l-x)m 2 -x(l-x)u' 
*2> *• *) = TTT^VtT^ I log 2 f ^ - x V log 2 f ^ - x) + 1 log 2 ^ " 



AV 2 (u,m 2 ,M 2 ) 1 * \2u J ^ \2u J 2 * 2AV 2 (u, m 2 ,M 2 ) 
2vr 2 /2A 1 / 2 ( u ,m 2 ,M 2 )\ ( 2xu 



+ + Lb [~ u:. 2 2,„ ' ) + Lb UaV^u, m\Ml) ) ) ifX< 5' 



7T 2 



, A 1 / 2 (u,m 2 ,M 2 ) u+-2xu ( ut - 2xu 

+ log i ^ log - Ll 2 ' 



u 2xu-^ 2 V2A 1 / 2 (u,m 2 ,M 2 ) 

;u — 

V2A 1 / 2 ( u ,m 2 ,M 2 ) 



. / 2xu -u l2 \ I u 12 ut 2 



*§W,*) = TTT^-VttJ log 2 fx - #V log 2 fx - ^ ) - I log 2 



A 1 / 2 ( u ,m 2 ,M 2 ) | b V 2^7 b V 2«/ 2 b 2A 1 / 2 (u,m 2 ,M 2 ) 
vr 2 _ / u+-2xn \ / 2A 1 / 2 (n,m 2 ,M 2 ) 



+ Y" Li2 UAV 2 Km 2 ) M 2 )J- Li2 l 2xu-nr 2 ) f ^^7' 



u 



(s, t, x) = 4xm 2 M 2 - M*(l + 2x) - (2(1 - x)m 2 - M 2 ) (u - m 2 ) 



+ A ccx ((1 + x)M 2 - (1 - x) (u - m 2 )) . (B.39) 
B.2.12 1 meson, 1 nucleon, 2 photons 

1 /• d d fe 1 

11 (S, ) - . ^ (2?r)d ^ _ fc 2) ( m 2 _ (jfe _ A )2) ( M 2 _ ( g + fc )2) ( m 2 _ ( p _ fc )2) 

d d A: 1 



i 7i (2vr) d (m 2 - (k - q) 2 ) (m 2 - (k - q') 2 ) (M 2 - k 2 ) (m 2 - (S - fe) 5 
~ Wn l m 2 t (l°g^2-- L ) 5n (^m^J-^^t), 
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II-, ( 4m 2 arccos 

(st)= dx dy^—^l 2my/s^(s,t,x,y) 

J J 167f2 I {-y 2 A(s,m 2 ,M2) -4x(l-x)(l-y) 2 m 2 t} 3/2 



+ s n (s,y) 



sl^s, t, x, y) {-y 2 X (s, m 2 , M 2 ) - 4x(l - x)(l - y) 2 m 2 t} J ' 
a?! (s, *, x, y) = a - (1 - y) (s - m 2 + M 2 ) + (1 - yf [M 2 - x(l - x)t) , 

oo 

Pn(s,t) =P U (so, t)+ J RP 1 ( So ,t) + _— ° y _ ^ ^_ ^ - J&fot), 
#(s-(m + M^) 2 ) / _ t x 

imp -^ = - ^(..^.W) ( los ^ " / • (R40) 

B.3 Tensor decomposition 
B.3.1 Meson integrals 

1 i J, (2ir) d M 2 - k 2 ' 



i r d d k w 



AM A 11 
= —Mt) + {{Mf - M 2 ) J^t) + Aj — A,) , J£(t) = — J«(*), 

m 1 f A(M 2 ,M 2 t) A- -A- 1 

12t\ 3 t 3y ' t v 1 3, \ 288vr 2 t 

a. (t + M 2 -M 2 Y ( t + M 2 -M 2 ) 

4 2) (*) = "^ + ^«W- (B.41) 
B.3. 2 1 meson, 1 nucleon 

i[S> i Ji (2ir) d [M 2 — k 2 ) (m 2 — (S — A;) 2 ) 1 K h 
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4 1} ( S ) = 1 (A, + (s - m 2 + Ml) /,(,)) , lW(s) = & -^-h(s). (B.42) 
B.3.3 2 mesons, 1 nucleon 

jvvm = I f d ^ fe ^ 

21 * Ji ( 27T ) d (M 2 - k 2 ) (m 2 - (fc - A) 2 ) (m 2 - (p - k) 2 ) 

= Q»l£ (1) (t) + + ^{h K) - J,- (m 2 ) + (M 2 - M 2 ) 4j(t)}, 



4 /i (*) = o M 4 (1) (*) + ^4i(*), 



* Ji ( 27r ) d (M 2 - k 2 ) (M 2 — (k — A) 2 ) (m 2 - (p - A:) 2 ) 
= sTJ^ 2 \t) + Q»Q v I l i {3 \t) + A"A"ig (4) (t) + i (A"Q" + Q M A") 4? 5) (t), 

4 (2) (0 = 4(23^ M + K) + V t ' ffift) + 2 (M 2 + M 2 - t) 

M 2 - M 2 1 
+ ^^(/,(m 2 )-/,(m 2 ))} 

i r x(m 2 ,m 2 t) .. 

= 4{ 7i M + /j M + t L l2 ^ + HM 2 + M 2 - t) I^\t) 

M 2 - M 2 1 
+ ^^(/,(m 2 )-/,(m 2 ))} 

2(1 + f(t, m)) (m 2 A (M 2 , M 2 , t) + tM 2 M 2 ) - (4m 2 - t) [m 2 - M 2 ) 2 
+ 128vr 2 m 2 (4m 2 - t) t ' 

4< 2 >(t) = -^— } {2/, (m 2 ) - (4M 2 - t) i«(t) + 2 (2M 2 - t) J*« '(*)} 

= -i{ 2 /, (m 2 ) - (4M 2 - t) + 2 (2M 2 - t) 4 (1) (t)} 

| (l + f(t,m))(M*-m 2 (4M 2 -t)) 
647r 2 m 2 (4m 2 — t) 



TLf 2 — M 2 

+ J 4 (m 2 ) + (m 2 ) + 2{d - 1) (M 2 + M 2 - t) ^(t) + ! 2 (I, (m 2 ) - /, (m 2 )) 
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= 8(4^y {* K) + 4 K) - 2 K) + M) + v t J ^ (t) 

fl/f 2 _ /Lf 2 ^1 

+ 6 (M 2 + M| - t) I^\t) + -L—^. (J, (m 2 ) - I 5 (m 2 )) | 

2(1 + f(t, m)) (m 2 A (ilf 2 , M 2 , t) + tM 2 M 2 ) - (4m 2 - t) [m 2 - M 2 ) 2 

128vr 2 m 2 (4m 2 - i) 2 t 

4 (3) (0 = 4(d _ 2) | 4ro 2_ t) {^ K) " 2 ( d ~ 2 )4 1} K) - (4M 2 - t) 

+ 2(d-l) (2M 2 -t)4 (1) (t)} 

= 8(4m 2_ t) { 2J i K) - 4 4 x) M - (4^ 2 - 1) im + 6 (2M 2 - 1) i£W(t) 

(1 + /(t, m)) (jg - m 2 (4M 2 - t)) 
647r 2 m 2 (4m 2 — t) 2 

lT\t) = " 2(d - 3)/, (m 2 ) + (d - 2) (if) (m 2 ) + (m 2 )) 



+ 2 (M 2 + M 2 - t) 4«(t) + (d - 1)^L_L 1 (J, ( m 2 ) _ /. ( m 2 )) 

1 



r ii(4) 



+ l{(d-l)\ {Ml Mf,t) + 4(d - 2)tM 2 } 

= 1{ - 27, (m 2 ) + 2 (if) (m 2 ) + J?> (m 2 )) + 2 (M 2 + M 2 - t) ^(f) 

f + M 2 — M 2 1 1 

+ '- [h (m 2 ) - Jj (m 2 )) + I {3A (M 2 ,M 2 ,t) + 8iM 2 } /£(t)J 

2(1 + f(t, m)) (m 2 A (m 2 , M 2 , tj + iM 2 M 2 ) - (4m 2 - t) [m 2 - M 2 ) 2 

128vr 2 m 2 (4m 2 - t) t 2 ' 

_(4M 2 -(d-l)t)4(t)| 

= 1{ - 21, (m 2 ) + 4/f) (m 2 ) + 2 (2M 2 - 1) I^\t) - {AM 2 - 3t) I«{t) 

_ (1 + f{t, m)) (M[ - m 2 {AM 2 - t)) 
647r 2 m 2 (4m 2 - t) t ' 

4 (5) W = Yt M " M) + 7 (* + m2 " M i) /2 ? 5) (*) = ^ (*)■ 

Note that /2i^(*)> ^2i^(*) an< ^ ^2i^(*) are symmetric in i and j, but that i^i^W an d ^2i^(*) are 
not. 



(B.43) 
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B.3.4 1 meson, 2 nucleons 



l 



I 



I (2vr) d [Mf - k 2 ) (m 2 - (p - A;) 2 ) (m 2 - (p' - /c) 2 ) 



4 1) (*) = t-^(^K)+mMw) 



4m 2 - t 

1 /• d d A: Wk v 



i Jj (2ir) d (M 2 - k 2 ) {m 2 - (p - k) 2 ) (m 2 - (p' - k) 2 ) 
= 9^lf\t) + Q^Q u lf\t) + A^A u lf\t), 

= Ml (r A{t) _ jf)( t j\ M ? W + " 4m2 + *) m )} 



2 V AW / 64vr 2 m 2 (4m 2 - t) 



m 



-{Af?(a7i«(t)-^(t))+J« (m 2 )} 



+ 



2 (4m 2 - t) 
M 2 {M 2 + (M 2 - 4m 2 + t) f(t, m)} 



647r 2 m 2 (4m 2 — t) 

= ± (*? (4% - r A(t) ) - 4" (m «) 1 + ^W + W-^' + Q/ft.")} 

2i\ *V A w AV V 4 v ; J 64vr 2 m 2 (4m 2 - t) t 

^ W (0 = /1 (2) (*) = lf ] (t) = IT Q) = 0- (B.44) 



[ J " i A (2vr) d (M 2 - A; 2 ) (m 2 - ( Pl - k) 2 ) (m 2 - (E — A;) 2 ) ^ + A i B W , 

= p^ + Sy*, A M = E M - pf, Q 2 = 2m 2 + 2s - M 2 , A 2 = M 2 , Q ■ A = s - m 2 , 

= — ( s " m2 " M l ) J * K) " ( s ~ ™ 2 + M l ) W 

2A (s,m 2 ,M?J L 
+ (( s _ m 2 ) 2 _ m 2 ( S - m 2 + 2M 2 )) 4( S )} ; 

j| {2) (s) = t — r-{ (3s + m 2 - M 2 ) /i(a) - (s + 3m 2 - M ? 2 ) h (m 2 ) 

2A(s,m 2 ,M 2 J l 

- (s - m 2 ) (s + 3m 2 - 2M, 2 - M 2 ) l|(s)}. (B.45) 

For the definition of pi, p2, and Mj, see Sect. 1572.41 In these conventions (in particular for Mj), the 
u-channel case follows by E —* A, s — > u. 
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B.3.5 1 meson, 3 nucleons 

k^ 



k 2 ) (m 2 — (p — k) 2 ) (m 2 — (S — k) 2 ) (m 2 — (p' — k) 2 ) 
= Q»I l £\s,t) + (A + 2^4 2) ( s ,t) + A^ 3 3) (a,i), 

P ™ M = W^J){ t(S ~ U)Tm " ( ^ ( S ~ m2 + M ") ~ (* + A ccx)Acex) 

+ (t (a - m 2 ) (a - u) - AtM 2 M 2 + M 2 (t + A ccx ) 2 ) P 13 {s, t) 
+ 2<f 13 (a) ((a - m 2 ) (t - A ccx ) + (t + A ccx ) (M 2 - A ccx )) } 

4(A( S ,^) + s t) {^-^- 2 ( g -- 2 + M -)™ 
+ (( s _ m 2) (a - u) - M 2 (4M 2 - t)) J} 3 ( S , t)], 



I 



7(1) 



(a, t) = m 1 {t(s - u) (Il(t) + (a - m 2 ) lUs, t)) 



13 4N 13 (s,t) 

- (2t (s-m 2 + Ml) - (t + A CCX )A CCX ) F*(s) 

+ 2Sj 3 (s) ((a - m 2 ) (t - A ccx ) + (/ + A ccx ) (M 2 - A ccx )) } 

A ccx =o, <5- 3 =o (a - u) (Jltt) + (a - m 2 ) Iljs, t)) -2{s-m 2 + Ml) Ig(s) 

4(X(s,m 2 ,Ml) + st) 

^ ){8,t) = ~ 4N 13 \s,t) {* (4m ' " im " * (2 ( S + m 2 - Ml) + A ccx ) /f (a) 
+ t((s-m 2 ) (4m 2 -t) - M 2 (s-u)) P n (s,t) 
+ 2<f 13 (a) ((a + m 2 - M 2 ) t - A ccx (2m 2 - t)) } 

+ ((a-m 2 ) (4m 2 -t) -M 2 (a- U ))/i 3 (a,t)}, 

t} = ~4A^M) {* (W " *) pA{t) ~ ' (2 (S + ^ " + Accx ) W 

+ t (s - m 2 ) (4m 2 - t) lj 3 (s,t) + 2(57 3 (a) ((a + m 2 - Ml) t - A ccx (2m 2 - t)) } 

A ccx =o, «57 3 =0 (4m 2 - t) {I\(t) + (a - m 2 ) 7? 3 (s, tj) - 2 (a + m 2 - M 2 ) 4"(a) 

4(A(a,m 2 ,M 2 ) + at) 

^ 3)(s ' *) = ii^T) { ( 2 ( s + m ' " M *) + Accx ) #W " ( 4m2 " *) ^ 

- ((4m 2 -t)(s- m 2 ) -M 2 (s-u + Ml)) P 13 (s, t)} 

+ Mjj) { ( 4m2 " *) ( 4M2 " 1 ~ 2Accx ) " (s " n)2 

- A ccx (2 (a + m 2 - M 2 ) + A ccx ) } Ac ° x= ^ 513=0 0, 

^ M = ii^j) { (2 (s + m2 ~ Ml) + Accx) ^ {s) ~ (W " t} m 
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(4m 2 -t) {s-m 2 )r i3 (s,t)} 



+ 



{ (4m 2 - t) (4M 2 - t - 2A CCX ) 



2 (s + m 2 -Mf) + A r 



=o, <57 3 =o 



Nis(s, t)=t(X (s, m 2 , Ml) + st) + tA cex (s + m 2 - M 2 ) + m 2 A. 



2 A 2 

cex> 



n 



2 I '// 





TZ7T 

1 B 



for 7T p — > 7r°n, 
otherwise. 



(B.46) 



The corresponding decomposition for the u-channel reads 



I^(u,t) 



d d k 



i Jj (2vr) d (M 2 - k 2 ) (m 2 - (p - k) 2 ) (m 2 - (A - A;) 2 ) (m 2 - (jjf - k) 2 ) 



Q^\u,t) - (A + 2qri^'(u,t) - A^(u,t), 



i(2), 



(B.47) 



where the basis functions as defined in (|B.46P may be used. The minus in front of I l 13 2 \u,t) is 
due to the replacement q <-> —q' to arrive at the u-channel, whereas the second minus is caused by 
our parametrization of the pion mass difference occurring in the kinematics of the charge exchange 
reaction: contracting (|B.46P with Q^, (A + 2q) fl , and A M , yields 



/ Q 2 Q-(A + 2q) 

Q ■ (A + 2q) (A + 2q) 2 -A, 
-A cex 



t 



iif\ s ,t) 

\4 3) M), 

4(t) + ( s - m 2 )ii3( S) t)-i^ (s)+ ^ 3 ( S) | , 



(B.48) 



where we have used that Mj = M % for p\ = p in the s-channel. Supposing that the sign of J^g (u, t) 
is not changed as compared to l[ 3 3 \s,t), we obtain 



/ Q 2 -Q-{A + 2q) 

-Q • (A + 2g) (A + 2<7) 2 A cc 
\ A cex i 

/ r%(u)-5\ 3 (u) + M 2 I{ 3 (u,t) 

P A (t) + {u-m 2 )l\ 3 (u,t) 



(u,t) 
WML 



r i(2) 

! 13 

r i(3) 



(B.49) 



V 



<Ji3(«) 



Changing the sign of l\ 3 3 \u,t) and multiplying the third equation by (—1) ensures that both systems 
of equations are identical up to s <-> u. 



B.3.6 1 meson, 1 photon 



1 



d d k 



k 11 ' 



i y : (27r) e! (m 2 - A: 2 ) (M 2 - (q - k) 2 ) 



(i) 



10 



2M 2 ' 



(B.50) 
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B.3.7 2 mesons, 1 photon 



V& ® ~ i I (27TY (m 2 - fc2) ( M 2 - (q - k) 2 ) (M 2 - (q> - k) 2 ) " (A + W V *®> 
B.3.8 1 meson, 1 nucleon, 1 photon 

&vS ) (s)+^vS\s), 

--A»v£\u)+p»vg\u), 
--A»v}t\u)+p'>*vg\u), 



d d k 








(27r) d (m 2 


— fc 2 ) (to 2 — 


(p-fc) 2 )(M 2 - 


(g + k) 2 ) 


d d A: 








(27r) d (to 2 


— k 2 ) (to 2 — 


(p> - A;) 2 ) (M 2 - 


-{q' + k) 2 ) 










(27r) d ( m 2 


— fe 2 ) (to 2 — 


(p-fc) 2 )(M 2 - 


(q> - k) 2 ) 


d d A; 









i J l {2ir) d (m 2 - k 2 ) (to 2 - (p 1 - k) 2 ) (M 2 - (q - k) 2 ) 



r (i) M = (s-m 2 - Ml) Vio - (s + m 2 - M 2 ) I n (s) 

X(s,m 2 ,Ml) 
, m _ 2sI n {s)-{s-m 2 + Ml)V w 



^= X(s,m 2 ,M 2 ) " ■ (B.52) 

S M j4 (l) (s)+p ^( 2 ) (s)) 

= AM5i ) (u)+^4 2 1 ) H, 



4?w 









1 






7i (2vr) d (m 2 - k 2 ) (to 2 - 


(S 


- A:) 2 ) (M 2 - 


(q-k) 2 ) 




/• d d A; 




1 




i , 


h (2vr) d (m 2 - A: 2 ) (m 2 - 


(S 


- fc) 2 ) (M 2 - 


(q> - k) 2 ) 




f d d k 




1 






h {2txY (m 2 - k 2 ) (to 2 - 


(A 


- k?) (M 2 - 


(q> + k) 2 ) 




r d d k 




1 






h {2txY (m 2 - k 2 ) (to 2 - 


(A 


- k) 2 ) (Ml - 


(q + kf) 




{s-m 2 - Ml) (Vio + (• 


s — 


m 2 ) An(s)) - 


- (s + m 2 



X(s,m 2 ,Ml) 



(2) (s - m 2 + Ml) (V 10 + (s- to 2 ) A n (s)) + (s + to 2 - M 2 ) (m 2 ) - 2s/ 7 (s) 

An (S) " A(s,to 2 ,M 2 ) • (B - 53) 

B.3.9 2 mesons, 1 nucleon, 1 photon 



/la- 



r ^ 1 f d ' 



(m 2 - k 2 ) (to 2 - (S - k) 2 ) (Ml - (q - k) 2 ) (Ml - (q> - k) 2 ) 
\s,t) + (^ + 2qYA { S(s,t), 
d k W 



21 ' ' i (m 2 - /c 2 ) (to 2 - (A - k) 2 ) (Ml - (q + k) 2 ) (Ml - (q' + k) 2 ) 
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Q"41w)-(A + 2?)> 4 4i (u,t), 



WW = 4(M..m',«|) + rt ) < 2 < S - m + - ( » - 



-(4Af;-*)((s-m 2 )A 21 (s,() + 

^ M = 4(M..m^) + rt ) { (4 '" 2 - ') - 2 (S + ^ - M -» 

+ (s - «) ((s - m 2 ) Aai(j,i) + Vioft)) }. (B.54) 



B.3.10 1 meson, 2 nucleons, 1 photon 

The tensor decomposition for the s-channel may be chosen as 

1 f d d k k» 



i Jj (2Tr) d (m^ - k 2 ) (m 2 - (p - A;) 2 ) (m 2 - (E - A:) 2 ) (M 2 - - A;) 2 ) 

= QMg( S ,i) + (A + 2g)Mg ) ( S ,t) + AMg^s.t), 
1 /" d d A; A^ 

7 / (2vr) d (m 2 - A; 2 ) (m 2 - (p' - k) 2 ) (m 2 - (S - A:) 2 ) (M 2 - (q - k) 2 ) 

= Q»A^(s,t) + (A + 2q)»A$(s,t) - AMg^a.t), (B.55) 

which for the u-channel becomes 

1 f d d k kf 1 



i J l (2vr) d (m 2 - A; 2 ) (m 2 - (p - A:) 2 ) (m 2 - (A - A:) 2 ) (M 2 - ( g + A;) 2 ) 
= Q"Ag («, i) - (A + 2qYA^ (u, t) + A"4? (u, t) , 
1 r d d k k» 

1 h {2ir) d (m 2 - A; 2 ) (m 2 - (p' - A:) 2 ) (m 2 - (A - A:) 2 ) (M 2 - + A:) 2 ) 

= QM$(u,i) - (A + 2q)>*A ( $(u,t) - A^? («,*)■ ( B - 56 ) 

Since only the first version of (IB.56j) contributes to the charge exchange reaction, A cex is incorporated 
according to the prescription in Sect. lB.2.111 The pion mass which has to be inserted into the definition 
of l]j(s) is Mj = \J M 2 — A ccx . In the present context, we thus deviate from the original definition 
in Sect. [BXil 

A " (fl,t) = 4iV 12 1 ( s ,t) { (t ( u ~ m2 + M ') + A - ( s " m " ~ M ')) # («) 
+ (t (s - m 2 + M 2 ) + A cex (u-m 2 - M 2 )) i^'to 

- (t (s - m 2 + M 2 ) + A cex (s - m 2 - M 2 )) A u («) 

- (t (u-m 2 + M 2 ) +A ccx (u-m 2 -M 2 )) ((s - m 2 ) A l2 (s,t) + V n {u)) } 



A cex =0 



.| ( u _ m 2 + m 2 ) I^(u) + {s-m 2 + Ml) Tg(a) 



A(X(s,m 2 ,M 2 )+st) 
- (s-m 2 + M 2 ) An(s)- (u-m 2 + M 2 ) ((s-m 2 ) Ai 2 (s, t) + Vn( 

^ 2 (S ' * } = 4iV 12 1 ( s ,t) { (t ( M + m ' " M ') + A - (* " 2 ™ 2 )) 7 * : («) 
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- (t (s + m 2 - Ml) + A cex (t - 2m 2 )) Pj (s) + (t (s + m 2 - M 2 ) + 2m 2 A cex ) An (a) 
-(t{u + m 2 -Ml)+ 2m 2 A ccx ) ( (a - m 2 ) A 12 (s,t) + Vu («)) } 

+ (a + m 2 - Ml) A n (s) - (u + m 2 - M 2 ) ( (a - m 2 ) A 12 (s,t) + V u (t 



= ~^{^b'(«) " i#00 " An( S ) + (a- m 2 ) A 12 (s,t) + Vh(u)} - ^A®(s,t) 

A = =0 " 4» " ^liW + (* - ™ 2 ) ^M) + Vii(«)}, 

N 12 (s, t)=t(X (s, m 2 , M 2 ) + st) + A cex ((a + m 2 — M 2 ) i + m 2 A cex ) . (B.57) 

B.3.11 1 meson, 1 nucleon, 2 photons 

P . (s t) = i r^L k i 

lU ' ; ijj {2n) d (m 2 - k 2 ) (m 2 — (k — A) 2 ) (M 2 - ( g + /c) 2 ) (m 2 - (p - A;) 2 ) 
= -<ffii(a,i) + Q^M) + (A + 2g)"P 1 ( 1 2) (a,t), 

j» ( f) = 1 f d d k W 

11 {Ul ' % h (2vr) d (m 2 - A: 2 ) (m 2 - (fc - A) 2 ) (M 2 - (g' - fc) 2 ) (m 2 - (p - fc) 2 ) 

= ^Pii(M) + Q^PgW) - (A + 2g)"P 1 ( 1 2) (M), 

^M) = - 4 (A(^ 2 lM 2 ) + S t) {( 5 " + ^ " *) 

_ ( s _ m 2 + m 2 ) (tP u (a,t) + 2V u (a)) }, 

^(s,t) = ^ + 4(A(g , m2 ; M , ) + st) {(^ 2 " + (- " «)flo(*) 

- (s + m 2 — M 2 ) (tP 11 (s,t) + 2V 11 (s))y (B.58) 

B.4 Bremsstrahlung 

The following representation of the integrals appearing in (|2.77p is valid for s > (m + M w ) 2 : 

7 p z |p| Pp - |p| 



, . 2 s + m — M 2 s + m 2 -M 2 + A 1 / 2 (s,m 2 ,M 2 ) 
5pp(s) = ^2AV2( Sjm 2 >M 2) lo § 2nVi ' 

2 a - m 2 + M 2 a - m 2 + M 2 + A 1 / 2 (a, m 2 , M 2 ) 
9qq[S) ~ M 2 AV 2 ( s , m 2 ,M 2 ) g 2M^7i ' 



5pg(«) = 2 



/■ dx s + (2x-l)(m 2 -M 2 ) 



Spg(x) ^ ((a + (2a; - l)(m 2 - M 2 )) 2 - 4 SSpg (x) 

s + (2a;-l)(m 2 -M 2 ) + J (s + (2s - l)(m 2 - M 2 )) 2 - 4ss pg (s) 
x log — . 

2y/SS pq (x) 
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9pp 



Sp q (x) = x(l - x)s + (2x - l)(xm 2 - (1 - x)M%) 
= -x(l -x)(u + t)+ xm 2 + (1 - x)M 2 , 

} dx s + m 2 -M 2 s + m 2 -M 2 + Js 2 ppf (x) 

(S,t)=2 / = ; r , log , 

J m *- x (l- x )t ^Jxj 2 W™ 2 - x(l - xjt 



o 

l 



r dx 3 _ m 2 + M 2 ^ s - m 2 + Ml + ^ s 2 qql {x 

g ^ t} - 2 J Af?-x(l-x)t ^) l0g 2^Ml-x(l-x)t 

s 2 p , = s 2 q , = X(s,m 2 ,M 2 ) +Astx(l -x), 

( } dx s + (2x-l)(m 2 -Ml) 

g pq >{s,t) = 2 / — — ' 

j s pq > [X) ( s + (2a; _ l) ( m 2 _ m%) ) - 4ss p9 , (x) 



s + (2x - l)(m 2 - Ml) + J(s + (2x - l)(m 2 - M 2 )) 2 - 4ss M /(x) 
x log — ^=^= , 

2y / SS p9 /(x) 

s p<j'( x ) = — x)u + xm 2 + (1 — x)Ml 

= x{\ - x)(s + t) + (2x - \){xm 2 - (1 - x)Ml). (B.59) 

C Contributions of individual diagrams 
C.l Strong diagrams 

(-i) 

KM) = 4M 2 ^±4^(«) " (3* + m 2 ) 
s — m z 

. . s 2 + 10sm 2 + 5m 4 2 s 2 + 6sm 2 + to 4 r m , , . 
BIM = {s _ m2f M 2 Us) 4 >(,), (C.l) 

= ^ (2^( S ) + <( S )) , 5 S V( S ) = -^j (2B* Sl {s) + <(,)) , 

< p (*,t) = g£ (2^(n) + <(«)) , B^(s,t) = (2B-(u) + *£(«)) , 
\/2mg A 



(«2) 



^f(M) = (2^(n) +<(n) - 24*00 - , 

*.T(M) = ^ (2^(«) + <(«) + 25- ( S ) + . 

A£ (s,t) = 4 (s + m 2 ) - (3s + m 2 ) if } (s) + 4m 2 /f } (m 2 ) 



+ 8m 2 



2M 2 I i i{s) - 4 (s + m 2 ) - 2 (s - m 2 ) J^' (2) ( S ) 



M 2 /i(t) + 4m 2 /^ (1) (t) - (s - u) lf\t) + 32m 4 (s - m 2 ) f£\s, t), 



,2\ A 1 ), 
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B%(s,t) = (s + 3m 2 ) Ii(s) + 4m 2 4 1} (m 2 ) - (s + m 2 ) 



r(i), 



+ 4m 2 



M 2 ^(t) - 2/f )(t) + 2M 2 /^( S ) - 2 ( a + 3m 2 ) /g«( S ) 



2 (s - m 2 ) /^ 2 )(,)] - 16m 4 (M?P 13 (s,t) - 2 (s - m 2 ) 4 2 >(M)) , 



m<7 



\/2mg 4 
32F 4 



1/4 ^;M)-^<>,i), 



(<-(,,t)+^ o ( S ,t)) 



\/2mg 4 



32F 4 



(A^(u,t)+A^\u,t)), 



Bs?(s,t) = (Bf 2 *(s,t) + Bf 2 *°(s,t)) + ^ (Bf 2 «{u,t)+Bf/{u,t)) . (C.2) 



A%(s) = (s- m 2 ) 4 1} ( S ) - 2A, - 8m 2 (m 2 /^( S ) - (s - m 2 ) /g (2) ( S )) , 
B%(s) = -M 2 / 4 ( S ) + (a - m 2 ) - 4m 2 l|«( S )) 

4 " 2 (A, + (« + 3m 2 ) (M 2 /g( S ) - ( a - m 2 ) jg< 2 >(*))) , 



+ 



s — m? 



16F 4 
16F 4 

^ ^BJ^u) + <» - <^°(n)) . (C.3) 



^ S4 ( S ) = -2M 2 /,( S ) + ( S -m 2 )/W( S ), 
Bl 4 (s) =-(s + m 2 ) /«(,) + Mf-^I t (s), 

< P {s,t) = ^ (A* («) +<(n)) , < P (M) = -^4 (W +<(u)) , 
*) = - ^ (a^ W + <(,)) + (3Al(u) + Af M ) , 
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^T(M) = (3Bl(s) + <(,)) - ^ (35- („) + <(«)) . (C.4) 

^ = ^ = ^ 5 +<, i? s v p ( S ) = J Br 5 ( S )+<( S ), 

= 0, £- x ( S ,t) = (3SJ 6 («) + <( S ) + 3Bl(u) + <(u)) 

(«e) 

A%(s) = (s- m 2 ) (-21, ( S ) + / + 8m 2 /^ (1) ( S )) , 
= —Mfli (m 2 ) - 2A, - Sm 2 M*I l £(s) 

+ 2 ( S - m 2 ) (j« (s) — (,) + 4m 2 (,) + > (,)) ) , 

< P (s, t) = ^ (A-( S ) + <») , <^ ( S , t) = ^ - <») , 

= ^ (- B sf(s) + - + ^(u)) . (C.6) 

(«7) 

= ^(--«)i5 S) (t), 
^(M) = ^ (A, - 4m 2 (/W ( ro 2) + M?P A {t) - 2lf\t))) , 
A*;r(8,t) = -A^(s,t) = Af 7 (s,t) - 2A* r (s,t), 
Bl p {s,t) = -B^(s,t) = B£(s,t) - 2B* r (s,t), 

A™(8,t) = y/2A? 7 (s,t), B^{s,t) = y/2Bf 7 {s,t). (C.7) 

A-K-p _ ATT + p _ A CCX _ r, 

Bl P (t) = 4l4 1 J(t), B^(t) = -±4 1 J(t), B^t) = -^jU(t). (C.8) 



= -B* 5 (u)-Bf 5 (u), 

(C.5) 
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(-9) 



<(-) = m[ \ F f ] I?{s), Bl g (s) = ^ (-4M 2 /,( S ) -A l + 4(s- m*) /«(,)) , 
AlP(s,t) = 2Af(s) + Al(s) + Al(u), BT'(a,t) = 2Bf(s)+Bl(s) - Bl(u), 



sg \"i "/ sg \ / ' sg\ / ' sg V / ' sg V ' / sg V / 1 sg \ / sg ' 

l£*(*,t) = A* (a) + 2<(u) + A* («), B^(s,t) = B£00 - 2<(u) - B£( 



A- x ( S ,t) = -V2K„( S ) - A* («)) , 5 s c 9 ex (M) = -V2(SJ 9 ( S ) + B* («)) . (C.9) 

(-10) 

ATT-p _ ATT+p _ A CCX -Q R i (l) ~ 92 A S + 3m2 

A s w - ^sio - ^ 4 ^r, ^sio - U ' ^siol-J - gp4^ s - m 2 ' 

B£» = -4Br M (5), < P M) = 45^ (u), 
B™(s,t) = V2 (BIJs) + Bfjs) + B£» + <») . (CIO) 



(-11) 



A iT-p _ Air+p _ rng 2 Ml , 2 , _ 
A sn - A Sll - pi j tt \m ) , A sn - u, 



= ~< P = (™ 2 ) , 5- = ^ ( M 2 /. (m 2 ) + M^o (m 2 )) . (C.ll) 

(-12) 



= (* - ^') {W^WW - J^(t)} + ^{±8m 2 (,-u) ^(t) 

- (t + 4 (M 2 - M 2 )) (4m 2 / 2 7 (1) (t) - J^t)) + 2M 2 / W (m 2 ) }, 

= -< P (i) = - Jjlto) + 4m 2 / 2 7 (2) (t)}, (C.12) 

= ^^(--n)/ 2 T° (3) (t), flS(t) = ^{*^(*) + W^w}. 



(-13) 



i i / i^ 1 /) _ ircx _ n P 7T i, _ pTT 1 / _ 4A7T + A^o x _ \/2 (3A,r + 2 A^o ) 

^SIS — ^513 _ ^13 — U > n S!3 — n S!3 — ^pi ' U S13 ~ g^i ' K^' 1 ^) 



(-14) 

= < p = A Z = = < p = B Z = o- (c.i4) 
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Vector-type electromagnetic diagrams 



= &A~ v ™ + 4 ( s - ™ 2 - M D v ^ s ) + 4 ( s - ™ 2 ) ^w}- 



acck rjccx n 

A V2 — n V2 ~ U - 



e 2 



K7 = K V = A IT = B vT = 0. ^a" P (*) = ~< P (t) =~ T2 (2m 2 - t) P A (t) 

An~p AK + P 4Cex rjcex r\ 

^1)4 ~~ ^4 ~~ ^4 ~~ n V4 — U ' 

£„V (t) = -Bf A V(t) = -J^{ (2M 2 - 1) (V 20 (t) - 2V$(t)) + V$ - V W }. 



ATT p _ ATT^p _ A CCX _ r>1T p _ f>TT^p _ E>CCX _ ri 



Al 6 p = At/ = AZ = 0, Bi; p = -Bfj = -2V2B™ = - jL ( 2 y 10 - 
Axial-type electromagnetic diagrams 



^:r p = < p = az = Bl/ = B£* = BZ = 0. 



A17 = At p = ~ 2 -^ (2V W - v£>) , A™ = 0, 
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(as) 



Ai 3 (s) = ~^f{(d - 2) (l 7 ( a ) - /«(,)) + 2(d - 4) ( a - m 2 - M 2 ) 7^( S ) 

- ((d - 4) ( a - m 2 ) + 4M 2 ) i£( s )} = ~^f{2 (1,(8) - #>(,)) - 4M 2 P B \s)} 
me 2 5 2 ( s _ m 2)2 me 2 g 2 M 2 (s - m 2 ) (s + m 2 - M 2 - 2sf u (s, m, Af,-)) 



32vr 2 F 2 s 2 87r 2 F 2 S A (s, m 2 , M 2 ) 

54 ( a ) = 2)/ 7 ( S ) -4m 2 /^'( S ) -4m 2 (d-4)/^' (1) ( S )} 



2eV 

"~p2- 



{/ 7 ( S )-2m 2 /^( S )} 



e 2 



+ 



5 2 (s - m 2 ) (s-m 2 - M 2 ) (s + m 2 - M 2 - 2s fn (s, m, Mjfj 
16vr 2 F 2 s\ (s,m 2 ,M 2 ) 



Al/(s) = Al 3 (s), BlP(s) = B- 3 (s), A^(s,t) = Al A (u), B^(s,t) = -B^(u), 



A^(s,t) = -^-A: 3 (u), B^{s,t) = ^-Bl z {u). (C.23) 



(04) 



o2„2 



O) = ip- {< M - 2 >M») - < d - 2)4"w} 

= («,(.) -/(■»(,)) + ^V J ((. + m »)»-*») 



o2„2 „2„2 , „ e 2^,2 ^2 _ m 4^ 



B a: P (s) = ^(d - 2) (/ 7 (*) - = ^ (W " 4^ 



64vr 2 F 2 s 2 ' 

A£*{8, t) = A£*(u), b£*{8, t) = -BZJ{u), Alf = KT = 0- (C24) 



(05) 



- _ ^ 2 g 2 
B ~ 



[V W - V$ +2 (s-m 2 - Ml) V n (s) + (s - m 2 ) (v^ (s) - vg\s)) }, 



= -^{2 (s + m 2 ) vff(s) - 4m 2 v£\s) 
^SfW = -^K?(s), B™(s) = -^B:;p(s), B^(s,t) = -BlP(u). (C.25) 



7T p 
15 
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(«e) 



A*~*{s,t) = ^^{Sm 2 (u - m 2 ) A ( ^(s,t) -2 (u-m 2 - Ml) V u (u) + h(u) 

- Ml (v£\u) + v£\u)) + (s-m 2 - 2Ml) (a u (s) - A$(s)) }, 

K, P {s,t) = ^r{v w - 8m 2 (u-m 2 - Ml) A 12 (s,t) + 16m 2 (u - m 2 ) A^(s,t) 

- 8m 2 r B 7T (s) - 2I n (u) + 2 (u-m 2 - Ml) V u (u) + 2M 2 V 1 ( 1 1) (u) 

- 2m 2 (2A u (s) - 3A$(s) - A<$(sj) }, 

A£*(a,t) = AlJ(u,t), B£'{a,t) = -B£*{u,t), A™(s,t) = -^AlJ(u,t), 
BlT(s,t) = ^BlJ(u,t) + 2 ^ e V ( W _ ^ (u) ) . (c.26) 



K?{s,t) = ^f{2tll(t) - I,(s) - d -^I?{s) - (AMI + (d~ 4) (s - m 2 )) rg{ t 

+ 2(d - 4) (s - m 2 - Ml) Il n{1) (s) + 4m 2 (s - m 2 ) lj 3 (s,t) 
+ 4m 2 (2 (u - m 2 ) + (s - m 2 ) (d - 4)) lj 3 {1) (s, t) } 

= ^f{2tP A (t) - I 7 (s) - lW{s) - AM 2 J^(s) + 4m 2 (s - m 2 ) P 13 (s,t) 
+ 8m 2 (u-m 2 )lj 3 {1) (s,t)} 

me 2 g 2 (s - m 2 ) ( s - m 2 (s - m 2 f - M$ - 4sM 2 /n (s,m,M n ) 



16tt 2 F 2 [ 4s 2 sA(s,m 2 ,M 2 ) 
2 (s - m 2 + Ml) f u (s, m, M n ) - (s - u) f(t, m) ' 



X(s,m 2 ,Ml) + st 
Bl P (s,t) = §J {2 (2m 2 - t) l\(t) + (d - 2) (/ 7 ( S ) - lU(s)] 

+ A(d - 4)m 2 (i2T(s) - 2/^ (1) (s)) + 8m 2 (/ - 2m 2 ) I^( 8 ,t) 

+ 4(6 - d)m 2 (Am 2 - t) lj 3 (1) {s,t) + A(d - 2)m 2 (2M 2 - t) lj 3 {2) (s, i)} 

= ^{2 (2m 2 - Jftf) + 2 (j 7 ( fl ) - /WO 



+ 8m 2 



(t - 2m 2 ) J^( fl ,t) + (4m 2 - t) l£\s,t) + (2M 2 - t) lj 3 {2 \s,t)] } 



eV (s - m 2 ) J s + m 2 m 2 (s - m 2 + M 2 ) - s (s - m 2 - Ml) f n (s, m, M w ) 



+ 



8vr 2 F 2 [ 8s 2 sA(s,m 2 ,M 2 ) 
2 (s + m 2 - Ml) f u (s, m, il^) - (4m 2 - t) f(t, m) 
2(A(s,m 2 ,M 2 ) + st) 



A£ p {*,t) = ATPfat), B^(s,t) = -B^(u,t), A™ = B™ = 0. (C.27) 
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(as) 



K~ P {s,t) = ~^f{2 (2M 2 - t) V 20 (t) - 2V 10 + J^(t) - 8m 2 (2M 2 - t) A$(s,t) 
+ 4m 2 (aW(s)+A<$(s)) -Am 2 I^\t)}, 

B£>{s,t) = (2M 2 - t) (v M (t) - 2vg\tj) - V 10 + + 4m 2 (A u (a) - A<i>(*)) 

- 4m 2 (2M 2 - i) (A 21 (s,t)- 2A$> (s, i)) } , 
<*>( S ,t) = < P (M), i?^( S ,t) = - J B a V(n,t), ^ = ^ = 0. (C.28) 



< P (*) = {-I,(s) + M 2 (2A u ( fl ) - aU(s)) } , 

i?:/( S ) = ^ {2F 10 - vff - 2m 2 (4V (,) + A<? {s)) } , 
< P («, t) = < P (^, = -Bl p (u), A c a f = B« = 0. (C.29) 



(aio) 



ATT V — p /ITT + P _ n 

^aio "aio aio aio ' 

^2_2„2 



A^(s,t) = V g 9 ( (u 2 - Qum 2 + 5m 4 + (d - 4) (u - m 2 ) 2 ) Uu) 

ai()V ' 2F 2 (u-m 2 ) 2 \ V v M ; ; 7V ; 

- 2 ^_% 2 { ^ - 6 — 2 + 5 - 4 ) 7 tH + ((« + m 2 ) 2 - 4u 2 ) #>(«)}, 
*) = A J? e2g \A (!6^ 4 " 2 (« - m 2 ) 2 + (d - 4) (4m 4 - (u + m 2 ) 2 )) /» 

+ w - 2) ((„ + + w) /<■)(„)} = ^;^-;' ) (o.3o) 



+ 



128vr 2 F 2 u 2 

e 5 



V2- 2 - 2 



2F 2 (u - m 2 



2 | (8m 4 - (n-m 2 ) 2 )/ 7 (n) + ((u + m 2 ) 2 + 4um 2 )4 1) (n)}. 



(an) 



ATT p D7T p /ITT+p OTT+p r, 

"an an an on u ' 

V2me 2 g 2 (d-2 



«M) = W + ( 2 (« " 3 - 2 + M2 °) + (d " 4) (« - ™ 2 )) i2T°(«) 

- (4 (n - m 2 ) + (d - 4) (2« + 2m 2 - M 2 ,)) /^° (1) (u) - dM 2 // (2) (u) - /»} 

= («) " W + 2 (« " 3m 2 + M 2 ) //(u) 
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- 4 (u - m 2 ) r/ (1) (u) - 4M2 // (2) (n)} 

V2me 2 g 2 (u - m 2 ) 2 f n + m 2 - M 2 - 2uf u Q, m, M^o) 

128tt 2 F 2 u 2 \ + " A(u,m 2 ,M 2 ) 

B^( a> t) = AF f { f_ 92 m2) { (2 (« " 3m 2 ) + (d - 4) (« - m 2 )) /» + 2(d - 2)m 2 Z«(u) 

+ (4m 2 (it - 5m 2 + 2M 2 ) + 4(d - 4)m 2 (u - m 2 )) I 77r ° (u) 

- 4m 2 (4 (u - m 2 ) + (d - 4) (3« + m 2 - M 2 ,)) /^° {1) (n) " 4dm 2 M 2 / 77r ° (2) (u)} 
= 2j? ^ e ^ 2) { (« - 3m 2 ) /» + 2m 2 / 7 1 )(n) + 2m 2 (u - 5m 2 + 2M 2 ) lf\u) 



- 8m 2 ((„ - m 2 ) «(„) + Ml^\u)) } - ^^ {^ 



2m 2 (3u + m 2 - M 2 ) - 2u (u + 3m 2 - M 2 ) f n (u, m, M^) 
+ A(.,m 2 ,M 2 ) >' 



V2 



= ™p 2 9 fco - V&> + 2/ 7 («) - (« - m 2 + 2M 2 ) A u («) + (u - m 2 ) A^(u)}, 

B ZA^) = ^f{ 2m2 - 4?H - 2A U («)) + 2y 10 - v« (c.32) 

(Vio - V$ ] + 2/ 7 (u) - 2M 2 A n (n)) |, < 

( ffl 13) 



Au 3 - ^aig - Aiia _ ^13 ~ u ' 



<;(M) = -±^{/ 7 ( U) + ^>M} 

m,p 2 n 2 (it. — m 2 ^ 2 

(C.33) 



_ _ ^meV / (1) x vgmeV (n-m 2 ) 2 

2jF 2 l^W +i 7 WJ+ 128vr 2 F 2 U 2 

= 2 , | (« " 3m 2 ) /» - (u + m 2 ) #>(„)) - ^^fe^ - 

2F 2 (u-m 2 )\ y ' 7V y v ; 7 W J 128vr 2 F 2 u 2 

C.4 Photon exchange 

A™ = B2, cx = for i = 1, . . . , 14. 

Pi Pi 11 

(Pi) 

^/ = < P = 0, i?;/(t) = -^(t) = -^J«(t). (C.34) 
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ATT p _ ATT^p _ D1T ?_ D1 r P _ f) 
A P2 - A P2 - °P2 - °P2 - U - 



A£*( 8 ,t) = -<*(*, t) = (s - u) / 2 7 (3) (t), 



= < p = 0, B£'(t) = -B£*(t) = ~J^(t). 



a;/ = = o, = = (m 2 ) 



A;;p(s, t) = A;^(s, t) = 8me 4 { (s - m 2 ) pg> ( s ,t) + {u-m 2 ) P$ (u, t) } , 
B;;v(s,t) = B;^(s,t) = -4e 4 {M 2 P n (s,t) - M 2 P n (u,t) 
- 2 (s - m 2 ) P^ ] (s,t) +2 (u- m 2 ) pfffat)}. 



at p = a* p — n 

^Pio pio u ' 

4e 
1 



BfJ{t) = -Bfj(t) = ^{{2Ml - t) (V 20 (t) - 2v£\t)) - V 10 + v£ ] } . 
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(Pll) 

<? = <? = o, = = x { 2Fl0 " • (c - 44) 

(P12) 

= = 0. b pV(*) = = x (2m2 " t} (c - 45) 

(Pis) 

^/(t) = A£"(t) = 4me 4 {/ 2 7(t) + [d- 2)I 2 7 (1) (t)} 

4 (4m 2 + tf( 
8ir 2 m (4m 2 — t) 

(pu) 



W{ig W + 2gW(,)}- ^+M , B -^<^0. (C.46) 



C.5 Bremsstrahlung 

Apart from m 2 terms in the denominators, the following expressions for the bremsstrahlung diagrams 
are still exact. is defined according to (|2.6ip . The amplitude for 7r + p — ► n + p follows from 

T; +p {p,q,p',q',k)=TfP(p,- q ',p',-q,k). (C.48) 

C.5.1 Vector-type 

(vi) 

T ,- P ,_ e i j(2q' + kr j'{2 q -kY \ _ y^e f (2q - fef 

- 2F 2 \ 2q'-k 2q-k J ' U1 ~ 2F 2 2q ■ k ' 1 J 

V9 n 7719 1 ' n„/ 1 ~T 



2F 2 \' 2p'-k 2p-k 

V2e f n 20' (p - j + m ) 7 M ' 
4F 2 \ 7 + 2p- fc 



^ = "T7^ 7^ + g ^ o ; z ' \- (C50) 



rpTT p fl 

V3 



2F 2 



T-icex fi 

V3 



V2e 

If 2 



■Y 



(C.51) 
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(p 2 ) 



(Pa) 

2eV 
(P - P') 
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